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ABSTRACT 

We have obtained Spitzer Space Telescope IRS 5.5 - 35 /im spectra of 59 main 
sequence stars that possess IRAS 60 /mi excess. The spectra of five objects 
possess spectral features that are well-modeled using micron-sized grains and 
silicates with crystalline mass fractions 0% - 80%, consistent with T-Tauri and 
Herbig AeBe stars. With the exception of r] Crv, these objects are young with 
ages <50 Myr. Our fits require the presence of a cool black body continuum, 
T gr = 80 - 200 K, in addition to hot, amorphous and crystalline silicates, T gr 
= 290 - 600 K, suggesting that multiple parent body belts are present in some 
debris disks, analogous to the asteroid and Kuiper belts in our solar system. The 
spectra for the majority of objects are featureless, suggesting that the emitting 
grains probably have radii a > 10 /mi. We have modeled the excess continua 
using a continuous disk with a uniform surface density distribution, expected 
if Poynting- Robert son and stellar wind drag are the dominant grain removal 
processes, and using a single temperature black body, expected if the dust is 
located in a narrow ring around the star. The IRS spectra of many objects are 
better modeled with a single temperature black body, suggesting that the disks 
possess inner holes. The distribution of grain temperatures, based on our black 
body fits, peaks at T gr = 110 - 120 K. Since the timescale for ice sublimation of 
micron-sized grains with T gr > 110 K is a fraction of a Myr, the lack of warmer 
material may be explained if the grains are icy. If planets dynamically clear the 
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central portions of debris disks, then the frequency of planets around other stars 
is probably high. We estimate that the majority of debris disk systems possess 
parent body masses, M PB < 1 M e . The low inferred parent body masses suggest 
that planet formation is an efficient process. 

Subject headings: stars: circumstellar matter — planetary systems: formation 

1. Introduction 

Giant planets are believed to form in circumstellar disks by either (1) growth of interstel- 
lar grains into planetary cores and subsequent accretion of gas into planetary atmospheres 
(Pollack et al. 1996) or (2) direct collapse via gravitational instabilities. In the core-accretion 
model, the first step toward building planets is the growth of small bodies from sub-micron 
interstellar grains into meter-sized bodies. Fits of the 10 jum silicate features observed toward 
pre-main sequence T-Tauri and Herbig AeBe stars suggest that the grains in these systems 
have grown to radii of a few /im at ages <few Myr (Forrest et al. 2004; Uchida et al. 2004; 
Honda et al. 2003; Bouwman et al. 2003). However, how and when micron-sized grains grow 
to kilometer-sized planetesimals, and the efficiency of this process are not well constrained. 
Mid-infrared spectroscopy of debris disks around main sequence stars may be used to detect 
larger grains and to infer the presence of kilometer-sized planetesimals at later ages. 

Obtaining high signal:noise ground-based mid-infrared spectroscopy of debris disks has 
been challenging because typical debris disks have 10 jum fluxes, Fj, < 1 Jy, making the 
majority of systems too faint to be studied spectroscopically from the ground and using 
ISO. The excellent sensitivity of the IRS (Houck et al. 2004) on the Spitzer Space Telescope 
(Werner et al. 2004) has enabled the spectroscopic study of large samples of debris disks, 
including all objects discovered with the IRAS and ISO satellites. We present the first 5 
- 40 /xm spectroscopic study of a large sample of debris disks to study the growth of large 
bodies in circumstellar disks and to elucidate the physical processes acting on micron-sized 
grains. 

Debris disks are dusty, gas-poor disks around main sequence stars (Backman & Paresce 
1993; Lagrange et al. 2000; Zuckerman 2001). Micron-sized dust grains are inferred to exist 
in these disks from measurements of their thermal emission at infrared through millimeter 
wavelengths. The estimated lifetimes for circumstellar dust grains due to radiation and cor- 
puscular stellar wind pressure (if the grains are small), sublimation (if the grains are icy), 
Poynting-Robertson and corpuscular stellar wind drag, and collisions are typically signifi- 
cantly smaller than the estimated ages for the stellar systems, suggesting that the grains 
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are replenished from a reservoir, such as sublimation of comets or collisions between parent 
bodies. 

The dust grains observed around debris disks may be produced when planets gravita- 
tionally perturb parent bodies producing collisions. In our solar system, the zodiacal dust 
possesses dust bands that thermally emit L IR / L* = 1CT 7 times the incident stellar light and 
that have orbital properties identical to asteroid families, suggesting that the dust bands 
are generated by collisions between asteroids in each family. Gravitational perturbations 
by Jupiter and other planets in our solar system are expected to cause the apsides and 
nodes of asteroid orbits to precess at different rates because of small differences in their or- 
bital parameters. This precession leads to asteroid collisions that generate the small grains 
observed in the dust bands. However, observed debris disks are typically 3-5 orders of 
magnitude more luminous than our zodical disk. Simulations of self-stirred disks suggest 
that the formation of icy planets, with radii 1000 - 3000 km, may trigger collisional cascades 
between the remaining nearby kilometer-sized planetesimals (Kenyon & Bromley 2004). Gi- 
ant planet migration, soon after formation, may also trigger collisions. The migration of the 
Jovian planets in our solar system during the Late Heavy Bombardment may have caused 
gravitational resonances to sweep through the main asteroid belt, sending asteroids into the 
inner solar system, producing the craters observed on old terrestrial planet surfaces (Strom 
et al. 2005). 

The presence of planets in debris disks may be inferred from the dynamical influence they 
exert on dust particles. Planets may produce central clearings in disks by gravitationally 
scattering dust grains out of the system, that are otherwise spiraling toward their orbit 
center under Poynting-Robertson and corpuscular stellar wind drag and by trapping grains 
into mean motion resonances (Liou & Zook 1999; Quillen & Thorndike 2002). Unfortunately, 
directly resolving all but the nearest debris disks at infrared and submillimeter wavelengths, 
to search for disk structure generated by planets, is challenging with current ground- and 
space-based telescopes. However, the radial distribution of dust in a debris disk may be 
inferred from modeling the infrared spectral energy distribution (SED), assuming that the 
system is azimuthally symmetric. 

Mid-infrared spectra of 7 nearby debris disks, obtained with Spitzer IRS, have revealed 
that the majority of these systems do not possess spectral features, suggesting that the grains 
probably have a > 10 /im (Jura et al. 2004; Sloan et al. 2004). Jura et al. (2004) and Sloan 
et al. (2004) model the excess continuua of these objects using a continuous disk with a 
uniform surface density, expected if Poynting-Robertson drag is the dominant dust removal 
mechanism, and using a single temperature grain model. The disks around two-thirds of 
the stars in their sample appear to possess inner truncations at 10 - 50 AU. One possible 
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explanation for the presence central clearings is a planet at the truncation distance that 
sweeps the inner regions of the disk clear; however, central clearings may also be produced 
by sublimation if the grains are icy or by radiation pressure if the disks are collisionally 
dominated and grains with sub-blow out sizes are removed by radiation pressure. The SED's 
of the remaining systems are consistent with a disk with a uniform surface density that may 
be produced when large grains, generated by collisions between parent bodies, spiral into 
their orbit center under the Poynting- Robertson effect (and stellar wind drag). 

We report the results of a Spitzer IRS study of 59 main sequence stars with published 
IRAS 60 fim excesses, building on initial results published by Jura et al. (2004). Jura et 
al. (2004) included LL spectra obtained for 19 debris disks. This study includes SL and LL 
or SH/LH data, depending on the IRAS 25 fim flux, and expands the target sample size. 
Refinements in the IRS flats have allowed us to improve the signal: noise at which the spectra 
are measured; therefore, we reanalyze objects in Jura et al. (2004) in addition to presenting 
new data. We list the targets for the full sample, along with their spectral types, distances, 
and published ages in Table 1. 



2. Sample Selection and Characterization 

The debris disks around Vega, Fomalhaut, e Eridani, and f3 Pictoris were initially dis- 
covered from the presence of strong IRAS 60 fxm and 100 /iin excesses, 10 - 100 times larger 
than expected from the photosphere alone (Backman & Paresce 1993). Studies comparing 
the IRAS fluxes with predictions for the photospheric emission of field stars subsequently 
discovered more than 100 debris disk candidates, corresponding to a disk fraction ~15% 
(Backman & Paresce 1993; Cote 1987; Mannings & Barlow 1998; Sadakane & Nishida 1986; 
Sylvester et al. 1996; Walker & Wolstencroft 1988). 

We selected 115 stars with IRAS 60 fim excesses, discovered in the studies listed above, 
for further study with the IRS. Approximately 25 of these targets are pre- main sequence 
Herbig Ae/Be stars that will be discussed in L. Keller et al. (2006, in preparation) and 
another 30 of these targets are extended in the IRS slit, suggesting that the excess emission 
is generated by interstellar grains, and will be described in B. Sargent et al. (2006, in 
preparation). The discovery of extended infrared excess around two objects (HR 1307, 
HR 2522) is consistent with coronagraphic imaging that revealed the presence of reflection 
nebulosities in these systems(Kalas & Graham 2002). Modeling of the scattered light images 
suggest that the dust grains are small, interstellar grains, located at distances 1,000 - 100,000 
AU, rather than large circumstellar grains at distances <100 AU from the central star (Kalas 
& Graham 2002). 
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We list the debris disks in our study with their stellar properties and reported IRAS 
25, 60, and 100 excesses in Table 1. The majority of the stars in our sample are nearby 
(within 150 pc), isolated main sequence stars; although a couple of objects may lie outside 
the local bubble and another nine objects are members of well-studied OB Associations or 
moving groups. For example, the star A Cas is a member of Cas-Tau (de Zeeuw et al. 1999) 
with an estimated age of 10 Myr (Bhatt 2000). The star HD 146897 is a member of Upper 
Scorpius and the stars HD 95086, G Cen, HD 110058, and HD 113766 are members of Lower 
Centaurus Crux in Sco-Cen (de Zeeuw et al. 1999) with estimated ages of 5 Myr and 16 Myr 
(Mamajek et al. 2002), respectively. The stars HR 7012, rj Tel, and HD 181327 are members 
of the (3 Pic moving group with estimated ages of ~12 Myr (Zuckerman et al. 2001). 

Since we would like to study the evolution of dust properties as a function of time, 
we estimate ages for as many stars in the sample as possible by fitting logg and T e ff to 
Schaller et al. (1992) isochrones for B-, A-type, and F-type stars. We estimate logg and 
T e ff (see Table 1) from mean General Catalogue of Photometic Data Stromgren photometry 
(Mermilliod et al. 1997) using the calibration of Napiwotzki et al. (1993) and the rotation 
correction of Figueras & Blasi (1998). B-, A-, and F-type stars which appear above the 
published isochrones are assigned ages of 1 Myr, 50 Myr, and 300 Myr, respectively, the ages 
for which the main sequence isochrones begin. We estimate uncertainties in our estimated 
ages of a factor of two. For comparison, we also list published isochronal and moving group 
ages in Table 1. For cases in which the moving group age is inconsistent with the isochronal 
age, we assume that the moving group ages are more accurate because stars above the zero 
age main sequence in the HR diagram may be either evolving down onto the main sequence 
or up and away from the main sequence. 



3. Observations 

We obtained IRS spectra of 59 main-sequence stars with previously reported IRAS 60 
fxm excesses (Backman & Paresce 1993; Cote 1987; Mannings & Barlow 1998; Sadakane & 
Nishida 1986; Sylvester et al. 1996; Walker & Wolstencroft 1988) with either (1) both the 
Short-Low (5.2 - 14.0 fim) and Long-Low (14.0 - 38.0 /im; A/AA ~ 90) modules or (2) 
the Short-Low, Short-High (9.9 - 19.6 /im), and Long-High (18.7 - 37.2 /im; A/AA ~ 600) 
modules. In order to avoid time-consuming peak-up on our relatively bright targets with 
accurately known positions, we operated the observatory in IRS spectral mapping mode 
where a 2 x 3 raster (spatial x dispersion) centered on the star is performed (Watson et al. 
2004). We carried out the bulk of the reduction and analysis of our spectra with the IRS 
team's SMART program (Higdon et al. 2004). 
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We estimate the stellar photospheric fluxes of our objects by minimum x 2 fitting pub- 
lished photometry from the literature to model stellar atmospheres, using only bandpasses 
with wavelengths shorter than 3 /im: TD 1 (Thompson et al. 1978), Johnson et al. (1966), 
and 2MASS (Cutri et al. 2003). For stars with spectral types earlier than K2V, we use 1993 
Kurucz stellar atmospheres; for stars later than K2V, we use Nextgen models. We assume 
that all stars have solar abundances and log g = 4.5 unless otherwise noted. We use the 
Cardelli et al. (1989) extinction law to estimate E(B — V) and list the extinctions (A v = 
3.1 E(B — V)) estimated from photosphere fitting in Table 1. We plot the SEDs for all of 
our objects in Figure 1. Ultra-violet through near- infrared photometry are shown with black 
symbols; IRS spectra are shown in red; MIPS photometry, where available (Bryden et al. 
2005; Rieke et al. 2005), are shown with blue error bars; IRAS photometry are shown with 
green error bars and upper limit symbols; and submillimeter photometry, where available 
(Wyatt et al. 2005, Greaves et al. 2004, Sheret et al. 2003, Holmes et al. 2003), are shown 
with black error bars and upper limit symbols. Our photosphere models are shown with a 
solid black line. 

We measure the flux of our objects in two photometric bands to search for excess 
emission from silicates (8.5-13 /im) and cold grains (30-34 /im). The calibration uncertainty 
in the fluxes is ~5% and the measured statistical uncertainties are listed in Table 2. We find 
25 sources without strong IRS excesses despite reported IRAS 60 /iin excesses (annotated 
in Table 2). The average fractional excess, (F v (measured)-^ (predicted)) /^(predicted) for 
these photospheric objects is slightly negative in the 8.5 - 13 band (-0.028±0.055) with 
a standard deviation consistent with IRS observations of nearby, solar-like stars (Beichman 
et al. 2006). The average fractional excess is slightly higher in the 30 - 34 /im band and 
possesses a larger standard deviation (0.08±0.13), indicating that some of these sources may 
possess weak excesses at the longest IRS wavelengths. For sources that are not saturated 
and do not possess excess in the 8.5 - 13 /im band, we also normalized the photosphere to 
the first 10 data points in the SL module to search more sensitively for excesses in the 30 
- 34 /im band. Values of this calibration factor deviated from unity on a star-by-star basis 
by less than 5%. Pinning the photosphere to the fluxes at the shortest wavelengths of the 
SL module produces average fractional excesses for the 25 stars without IRS excesses of 
-0.032±0.048 and 0.06±0.13 at 8.5-13 and 30-34 /im, respectively. In our analysis, we use 
photosphere models that are scaled to the first ten points of SL unless the source is saturated 
or possesses a 8.5 - 13 /im excess. 

The significantly larger IRAS beam may contain multiple sources and therefore be 
responsible for some of the discprepancy between the IRAS and Spitzer IRS results. For 
reference, the IRAS beam is 1' x 5' at 12 /im and 25 /im, 2' x 5' at 60 /im, and 4' x 5' at 100 
/im. By contrast, point sources observed with Spitzer MIPS at have a FWHM of ~6" at 24 
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/im and a FWHM of ~20" at 70 /im. The color-corrected IRAS 25 /im fluxes are significantly 
higher than the IRS 25 /iin fluxes for three objects (HR 2124, HR 6297, and HD 200800). 
HR 6297 has been observed with MIPS at 24 /xm. In this case, the discrepancy between 
the HR 6297 IRAS and IRS fluxes can definitely be explained by the presence of a second 
source in the IRAS beam. The IRAS Point Source Catalog 25 jum flux for HR 6297 is (not 
color-corrected) ~410 mJy (with a central position 0.4' away from the star) while the IRS 
25 /im flux is ~90 mJy. MIPS 24 jum images reveal two sources near HR 6297; one at the 
position of HR 6297 with F„(23.68 /im) = 82 mJy, consistent with the IRS spectrum, and 
another, brighter source 0.92' west of HR 6297 with F„(23.68 /im) = 207 mJy. The IRAS 
PSC flux is the sum of the two MIPS sources and is therefore an overestimate of the HR 
6297 25 /im flux. Unfortunately, the remaining two sources with possibly confused IRAS 25 
/im fluxes, have not been observed with MIPS. We searched the 2MASS catalog around the 
positions of HR 2124 and HD 200800 for additional sources of confusion. For comparison, 
we examined the 2MASS colors for the unexpected 25 /im source (located at 17:00:00.07 
-54:35:38.6 (J2000)) in the HR 6297 IRAS beam (the spurious source possesses J=6.69, H— 
5.73, and If =5. 38). We found a bright source with similarly red 2MASS colors, 0.9' away 
from HD 200800 (located at 21:09:13.83 -65:47:38.4 (J2000)) with J=6.40, H= 5.47, and 
K—5.17. However, we were not able to find an obvious candidate for the extra emission in 
the IRAS HR 2124 beam. There are 7 2MASS sources within 2' of HR 2124 with J - K > 
1 but all have K ~ 14-15, suggesting that they are probably too faint to be detected at 24 
/im despite their red colors. 

The presence of IRAS 60 /im excess does not necessarily imply the presence of 5 - 35 /im 
excess. Twenty-five objects in our sample (~42%) do not possess strong infrared excess at 5 
- 35 /im despite the identification of IRAS excess in the literature. The discrepancy between 
published IRAS results and our results may (1) reflect difficulties in either the search criteria 
or the data reduction in the original searches, (2) be due to the presence of cold dust that is 
not detected at IRS wavelengths, or (3) be due to source confusion in the IRAS beam. Eleven 
of the objects in our study possess IRAS 12 /im and 25 /im fluxes that are consistent with 
our photosphere models and IRAS upper limits or photosphere detections in the IRAS Point 
Source or Faint Source Catalogs at A > 60 /im (107 Psc, CC Eri, r 1 Eri, a For, HR 1686, 
5 Dor, HR 3220, FI Vir, HR 3862, G Cen, and /i Ara) despite reports of measured excess 
in the literature. The discprepancy in these cases may be due to difficulties in the search 
criteria or data reduction. Our sample is drawn from a several IRAS searches for debris 
disks that used different criteria to establish whether a source possessed an infrared excess. 
Cote (1986) selected sources with V-[12], V-[25],V-[60] >0.5; Sadakane & Nishida (1988) 
selected sources with [12]-[60] > 1; and Walker & Wolstencroft (1988) selected sources with 
F^(60 /im)/F„(100 /im) = 0.8-2.0, similar to that observed toward (3 Pic, Vega, Fomalhaut, 
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and e Eridani. 

The dust in many of the undetected systems may be cold. The IRS spectra for five 
objects are photospheric over the bulk of the IRS wavelength range but may be rising at 
A ~35 /im although with insufficient signahnoise to model these systems in detail (k Lep, 
HR 2124, XZ Lep, 7 Tra, and HD 200800). The dust in some systems may be so cold 
that excess emission can not be easily detected above the photosphere at IRS wavelengths. 
Ten sources possess IRAS 60 /im flux excesses but no evidence for significant IRS excess 
at A < 35 pm (r Cet, 7 Dor, 32 Ori, ^ 5 Aur, 5 Vel, HR 5236, a Boo, 41 Ser, a Her, 
HD 221354). If the dust in these systems has T^ust < 60 K, then the infrared excess at 
the longest IRS wavelengths could be small and changing slowly compared to the stellar 
photosphere, making the excess difficult to infer. For r Cet, our measured and predicted 
fluxes, F„(8.5-13 /im) = 7470 and 7710 mJy, respectively, and i^(30-34 /im) = 786 and 829 
mJy, respectively, suggesting that this source possesses -3% and -5% excesses at 8.5-13 and 
30-34 /im, respectively. However, the 60 /im - 850 /im excesses observed toward r Cet are 
easier to detect because the photosphere is negligible at these wavelengths; these excesses 
are fit by a modified 60 K black body (Greaves et al. 2004), implying that the slope of the 20 
/tm - 35 /im spectrum should be ~4% higher than without the excess. Alternately, the IRAS 
flux for some objects may include extragalactic or other stellar sources that were included 
in the large IRAS beam. 



4. Grain Composition 

Five objects in our sample possess 10 /im and/or 20 /im spectral features: HR 3927, 
r) Crv, HD 113766, HR 7012, and rj Tel. We model the excess emission for these objects 
assuming that the observed features are generated by amorphous olivine and pyroxene, and 
crystalline forsterite, enstatite, and silica. We model the remaining continuum emission 
using amorphous carbon and one (or two) black body distributions. We do not observe PAH 
emission features toward any of the objects in our sample. 

A lower limit to the size of dust grains orbiting a star can be found by balancing the 
force due to radiation pressure with the force due to gravity. For small grains with radius a, 
the force due to radiation pressure overcomes gravity for solid particles larger than 

dmin,o = SL+Qpr/fi&KGMtCPs) (1) 

(Artymowicz 1988) where and M* are the stellar luminosity and mass, Q pr is the radiation 
pressure coupling coefficient, and p s is the density of an individual grain. We assume that 
the bulk density of silicate, amorphous carbon, and silica grains are p s = 3.3, 2.5, and 2.3 g 
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cm~ 3 , respectively. Since radiation from A-type (HR 3927, HR 7012, and rj Tel) and F-type 
(r) Crv and HD 113766) stars is dominated by optical and ultraviolet light, we expect that 
27ra/A ^> 1 and therefore the effective cross section of the grains can be approximated by 
their geometric cross section so Q pr ~ 1. We estimate the stellar mass by fitting our inferred 
T e ff and log g to the Schaller et al. (1992) isochrones for A-type stars and by fitting our 
inferred T e ff and L* to Siess et al. (2000) isochrones for F-type stars. We list the minimum 
grain sizes, a minj0 , for solid silicate, amorphous carbon, and silica grains around each object 
if the species is required in the fit in Table 3. If the grains are porous (with a volume 
fraction, / = V vac /V to t > 0), then the minimum-sized grains estimated in equation (1) must 
be modified to account for the vacuum fraction: 

a m inU) = yztj ( 2 ) 

In fitting the 5 - 40 /xm spectra of HR 3927, r] Crv, HD 113766, HR 7012, and rj Tel, we 
assume that the grains are spheres with radius, a > a m i n . We list the solid angles subtended 
by each dust population Q, their mass m and temperature, the radii of each grain population, 
and its vacuum volume fraction in Table 3. We infer absorption coefficients, Q a bsW, using 
optical constants published in the literature for amorphous olivine (MgFeSi04, Dorschner et 
al. 1995), amorphous pyroxene (Mgo.5Feo.5Si03 and Mgo.sFeo^SiOs, Dorschner et al. 1995), 
crystalline forsterite (Mg L9 Feo.iSi04, Fabian et al. 2001), crystalline enstatite (MgSi03, 
Jaeger et al. 1998), and amorphous carbon (Zubko et al. 1996), and Bruggeman Effective 
Medium Theory (Bohren & Huffman 1983). (Please see Sargent et al. 2006 for a more 
detailed description of QabsW estimates.) Non-zero vacuum volume fractions shift the peak 
positions of silicate features to longer wavelengths and broaden the features (Kessler-Silacci 
et al. 2006). 

Our fits to the HD 113766 and HR 7012 spectra require the presence of small sub-//m- 
sized grains (crystalline enstatite, silica, and amorphous carbon). In these cases, we modeled 
the sub-micron dust components using laboratory measured opacities of crushed forsterite 
(Mg 2 Si0 4 , Koike et al. 2003) and enstatite (Mg .7Fe . 3 SiO4, Chihara et al. 2002) and optical 
constants for silica (cristobalite, Simon & McMahon 1953) assuming that the distribution of 
grain shapes is well-described by the continuous distribution of ellipsoids (CDE; Fabian et al. 
2001). Both of these systems are young with estimated ages ~16 and ~12 Myr, respectively. 
The disks around around both stars are optically thin and possess high fractional infrared 
luminosities, L/r/L* = 0.015 and 10~ 3 , respectively. Since the estimate ages of these systems 
are ~16 and ~12 Myr, based on their membership in Lower Centaurus Crux in Sco-Cen and 
the (3 Pictoris Moving Group, these systems probably no longer contain bulk molecular gas. 
Therefore, we hypothesize that the small particles in each system may have been generated 
in a recent collisions between parent bodies. 
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We are confident in the identification of silicate species around HD 113766, HR 7012, 
and r) Crv, however, the exact dust mass in amorphous olivine, pyroxene and carbon is 
somewhat uncertain. For example, for HD 113766, the presence of sharply peaked spectral 
features at 10.0, 11.1, 12.0, 16, 19, and 23.5 are used to identify forsterite. For HR 
7012, the presence of spectral features at 9.3 and 10.5 fim are used to identify crystalline 
pyroxene. Once we fit the sharply peaked crystalline features, then we add amorphous olivine 
and pyroxene to fit the overall silicate feature; finally, we add amorphous carbon and one (or 
two) black bodies to fit the remaining continuum. The fit to the non-crystalline component 
of the spectrum is degenerate and may be fit with different ratios of amorphous silicates. 
In addition, the emissivity of amorphous carbon is approximately constant at 5 - 35 /xm, 
suggesting that a black body could be substituted for this component. In the case of HR 
7012, we use cristobalite because it produces a better fit to the 8.7 /xm shoulder of the 10 
/im silicate feature than alpha quartz; however, other silicates such as opal (hydrous silicate) 
may also provide good fits. 

The majority of debris disks with spectral features possess crystalline silicate emission 
features; only 77 Tel in our sample does not. We estimate crystalline silicate mass fractions 
of 76% and 0% for HR 7012 and 77 Tel (two ~12 Myr old members of the /3 Pic moving 
group), 4.1% for HD 113766 (a -16 Myr old member of Sco-Cen), 38% for HR 3927 (a -50 
Myr old field A- type main sequence star), and 31% for 77 Crv (a ~1 Gyr old field F-type 
main sequence star). The crystalline silicate fraction for HR 7012 appears extremely high; 
however, those of the other debris disks are consistent with measurements toward pre-main 
sequence T-Tauri and Herbig AeBe stars (Sargent et al. 2006; Bouwman et al. 2001). 
Therefore, we find no correlation between age and crystallinity for circumstellar silicates. 
The most dramatic example of the lack of correlation between crystalline silicate fraction 
and age is the disparity in the crystallinity of the dust around HR 7012 and 77 Tel, both 
A-type main sequence members of the (3 Pic moving group (Zuckerman et al. 2001). 

Our solar system is believed to possess two populations of dust grains produced by 
collisions in two distinctive small body belts: (1) the zodiacal dust located at D = 2 - 4 
AU with T gr = 270 K (Reach et al. 2003), generated by collisions between asteroids and (2) 
another population of cooler dust at D = 30 - 50 AU, with T gr = 50 - 60 K, that is believed 
to be generated by collisions between objects in the Kuiper Belt; this population of dust has 
not been detected directly thus far. Each of our systems with spectral features possesses 
at least two distinctive dust populations: a hot population with T gr = 290 - 600 K that is 
inferred to exist from detailed fitting of amorphous and crystalline silicate features, and a 
cooler population with T gr = 80 - 200 K that is inferred to exist from black body fits to the 
excess continuum. In addition, 77 Crv possesses a third, even colder dust component with 
T gr = 40±5 K, inferred from SED fits to IRAS 12 - 100 //m and JCMT SCUBA 450 and 
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850 /zm photometry (Wyatt et al. 2005). Our results establish the presence of multiple belts 
of small bodies in debris disks; previously, the evidence was much weaker. The "hot" dust 
component in these systems is only 2.5 - 3.6 x hotter than the "cold" dust component, while 
zodiacal dust is estimated to be 4 - 5x hotter than dust in the Kuiper Belt. HR 3927, rj Crv, 
HD 113766, HR 7012, and 77 Tel may possess asteroid and Kuiper belts in analogy with our 
solar system at larger distances from their central stars because of their higher luminosity. 
For rj Crv, the detection of three planetesimal belts with black body distances of 1.3, 11, and 
102 AU respectively may imply the presence of planets at 2.6, 22, and 204 AU (comparable 
to the distances of Mars, Saturn, and Sedna) if these belts are stirred by planets in the same 
way that the main asteroid belt is stirred by Jupiter in our solar system. Detailed dynamical 
models of this system are required to determine whether r\ Crv possesses multiple planets. 



5. Dust Properties 

The majority of infrared excess sources in our sample apparently lack spectral features, 
suggesting that the grains are too cold or too large (a > 10 /iin) to produce features. We fit, 
using x 2 minimization, the 5-35 /iin photosphere subtracted spectra for stars without silicate 
emission features, assuming that the grains are black bodies. The photosphere subtracted 
IRS spectra for all excess sources without spectral features are shown in Figure 3. The 
majority of our objects were observed in spectral mapping mode in which the slit is moved 
across the source, producing six independent spectra for each object. The error bars in this 
figure represent the difference between the two spectra in which the source is best centered. 
We overlaid our best fit single temperature black body model in blue. We list the best 
fitting grain temperatures, T gr , and fractional infrared luminosities, Li R jL if , inferred for 
each system assuming L IR = AilaTg r d 2 , where f2 is the solid angle subtended by the grains 
(in steradian) and d is the distance from the Sun to the central star, in Table 4 with the 
reduced x 2 f° r eac h fit- Since no significant excess is detected toward r Cet, despite a strong 
submillimeter excess characterized by T gr = 60 K, populations of cold dust may not be 
detected using IRS. Similarly, the temperatures inferred for cool dust populations using IRS 
spectra may be inaccurate. 

The minimum grain distance can be constrained from the grain temperature, T gr , as- 
suming that the dust particles act like black bodies. Black bodies in radiative equilibrium 
with a stellar source are located a distance 

D = k^) 2 R* (3) 
z i gr 

from the central star (Jura et al. 1998), where T e ff and are the effective temperature 
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of the stellar photosphere and the stellar radius. We estimate the stellar temperatures and 
absolute V-band magnitudes from Stromgren photometry using the calibration of Napiwotzki 
et al. (1993) corrected for rotation (Figueras & Blasi 1998) and infer luminosities using the 
bolometric correction from Flower (1996). We estimate R* assuming = AnRlaT^. The 
black body grain distances (listed in Table 4) range between 4 AU for the K1V star HD 
53143 and 72 AU for the B8Vn star A Cas. 

We estimate the average grain size of orbiting dust grains assuming that radiation 
pressure removes grains with radii < a miri)0 (please see equation 1) and that the dust grain 
size distribution is determined by collisional equilibrium (Greenberg & Nolan 1989) 



with p ~ 3.5, similar to that inferred for the interstellar medium based on interstellar extinc- 
tion curves, even though the grain-size distributions in these systems are not well constrained. 
Interplanetary and lunar size distributions with p sma ii = 2.2 - 3.7 for grains with a = 0.3 - 
2 jum and pi arg e = 2.0 for grains with a = 2 - 20 /im have been used to reproduce ISO 5 
- 16 /tm observations of zodiacal dust (Reach et al. 2003). If we weight by the number of 
particles, we estimate average grain sizes, < a > = 5/3 a mirij0 (see Table 4), assuming that 
the density of an individual grain, p s (= 2.5 g cm" 3 ). The majority of the stars in our sample 
are B-, A-, and F-type stars. (HD 53134 is the only object with infrared excess discussed 
here with spectral type later than F.) Since the radiation from these objects is dominated 
by optical and ultraviolet light, we expect that 2na/\ 3> 1 and therefore the effective cross 
section of the grains can be approximated by their geometric cross section so Q pr ~ 1. We 
estimate the stellar mass by fitting our inferred T e // and log g to the Schaller et al. (1992) 
isochrones for B- and A-type stars and by fitting our inferred T e // and to Siess et al. 
(2000) isochrones for the remaining stars. We find average grain sizes between 0.2 /im for 
HD 53143 and 25 /im, for the B8Vn star A Cas (see Table 4), similar to those inferred for 
the zodiacal dust in our solar system (Reach et al. 2003). 

We can estimate the minimum mass of dust around objects in our sample assuming that 
the particles have radius a minj0 ; if the grains are larger, then our estimate is a lower bound. 
If we assume a thin shell of dust at distance, D, from the star and if the particles are spheres 
of radius, a, and if the absorption cross section of the particles equals their geometric cross 
section, then the mass of dust is 



(Jura et al. 1995) where L IR is the luminosity of the dust. We estimate dust masses in 
micron-sized infrared emitting grains between 3.7xl0~ 8 M e for HD 53143 and 2.3xl0~ 4 



n(a)da = n Q a p da 



(4) 




p s D 2 a, 



•min,o 



(5) 
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M e for the A3IV/V star HR 1082 (see Table 4). If the grains possess a size distribution 
n G a~ 3 ' 5 with a maximum radius a max = 10 cm, as inferred from submillimeter observations 
(Zuckerman et al. 1995), we can estimate a the dust mass in larger grains. If the measured 
excess flux is F u (excess) and the black body flux for the excess is B v (excess), at frequency, 
u, then 

4 2 F u (excess) 

MlOcm = TTPv Q"min,oQ"max " ,-> / \ (P) 

3 v B u (excess) 
where d is the distance to our sun. We use the solid angle subtended by the dust grains, 
inferred from the minimum \ 2 black body fits to the excess, Q, to determine the flux to 
black body ratio, ^ v ^ xcess \ — Q anc l li s t the estimated larger grain dust masses in Table 4. 

J ' B v (excess) & & 



6. Dust Removal Mechanisms 



Dust grains in debris disks may be removed by a variety of processes such as radiation 
and corpuscular stellar wind pressure, ice sublimation, and collisions. Collisions may shatter 
parent bodies into small grains that are radiatively drive grains from the system. Larger 
grains in these high density environments may continue to collide until they reach sizes 
below the blow-out limit and are radiatively ejected. In lower density disks around B- 
and A-type main sequence stars, large grains may be subject to the Poynting-Robertson 
effect and may spiral in toward the central star. In the absence of planets in our solar 
system, Poynting-Robertson and corpuscular solar wind drag would determine the spatial 
distribution of dust. The discovery of debris disks around lower-mass solar-like and M-type 
stars has led to speculation that corpuscular stellar winds may contribute to grain removal 
in a manner analogous to radiation pressure and the Poynting-Robertson effect (Plavchan 
et al. 2005): (1) An outflowing corpuscular stellar wind produces a pressure on dust grains 
which overcomes the force due to gravity for small grains. (2) Particles orbiting the star are 
subject to a drag force produced when dust grains collide with atoms in the stellar wind. 
These collisions decrease the velocity of orbiting dust grains and therefore their angular 
momentum, causing them to spiral into the central star. 

The dominant grain removal process within a disk is dependent not only on the lumi- 
nosity of the central star but also on grain distance from the central star. For Fomalhaut, 
Backman & Paresce (1993) estimate that at 67 AU collisions to sizes below the blow-out 
limit (grains below the blow-out limit are quickly radiatively ejected from the system) are 
the most effective grain removal mechanism while at 1000 AU the Poynting-Robertson effect 
is the most efficient grain removal process. In Figure 4, we plot the sublimation lifetime, the 
Poynting-Robertson (and corpuscular stellar wind) drag lifetime, and the collision lifetime 
for average-sized grains around typical B5V, A5V, and F5V stars. Sublimation may quickly 



-14- 

remove icy grains in the innermost portions of the disk. At larger radii, collisions dominate 
grain destruction, and at the largest radii, where the disk has the lowest density, Poynting- 
Robertson and corpuscular stellar wind drag may dominate grain destruction. We estimate 
the lifetimes of average-sized grains and the parent body masses around stars in our sample 
if ice sublimation, Poynting-Roberton and corpuscular stellar wind drag, and collisions are 
each the dominant grain removal mechanism in the absence of other processes. 

6.1. Ice Sublimation 

One possibile explanation for the presence of central clearings, inferred from black body 
fits to the IRS spectra, is that the grains are icy and sublimate when they come too close 
to the central star. Although the dust grain composition can not be determined directly 
from spectral features, it may be inferred from the statistical grain properties in our sample. 
We plot the distribution of inferred grain temperatures in Figure 5. The estimated grain 
temperatures appear to cluster between 110 K and 130 K. Laboratory studies find that 
thermal desorption of water ice (H 2 0) from H 2 layers begins at temperatures of 120 K and 
is completed by 170 K (Fraser et al. 2001). Therefore, the peak in grain temperature at 110 
K to 120 K may suggest that the grains are icy and are beginning to sublimate. 

If the grains are icy, then sublimation may also remove grains from the disk. If the grain 
temperature is constant while the star remains on the main sequence, then we may write 
the following expression for the sublimation lifetime of an average grain. 

< a > p,iTg/ 2 e T ^ b ' /T ^ 
tsubl : 

(Jura et al. 1998) where d is the mass rate per surface area (= 3.8 x 10 8 g cm~ 2 s _1 K 1 / 2 , 
T su bi — 5530 K; Ford & Neufeld 2001) and pi = 1.5 g cm -3 if the grains are mostly ice 
with some refractory material mixed in, as expected for Kuiper Belt objetcs. The average 
grain radius, < a >, listed in Table 4 is computed assuming that the grains are composed 
of silicates with p s = 2.5 g cm -3 ; therefore, if the grains are mainly composed of ices, then 
the estimated average grain radii are an underestimate and the values in this table should 
be multiplied by 5/3. We list the sublimation lifetimes for grains in Table 5 assuming that 
they are icy. The sublimation lifetimes are sensitively dependent on grain temperature. For 
example, 3.5 pm grains around HR 1082 with an estimated T gr = 70 K, have a sublimation 
lifetime, T su m = 1.3 x 10 7 Gyr while 16 pm grains around HR 6211 with an estimated T gr 
= 160 K, have a sublimation lifetime, T su ^ = 7.4 minutes. In systems with T gr > 100 K, the 
sublimation lifetime is the shortest lifetime by more than an order of magnitude. 



(7) 
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6.2. Poynting-Robertson Drag 



If Poynting-Robertson (PR) drag is the dominant grain removal mechanism, then grains 
spiral in from the radii at which they are created toward their orbit center, creating a 
continuous disk with uniform surface density and a 1/fl volume density. In this model, 
the inner radius of the disk coincides with the stellar radius unless the grains sublimate, or 
are dynamically ejected by a massive body interior to the radius at which the grains are 
produced. For an optically thin, gas-free disk whose particle density, n, varies as D~ q , the 
infrared spectrum should be well-described by the function 



(Jura et al. 1998) where the absorption coefficient for the grains, Q a bs v p and K 1 and 
K2 are constants. The first term describes the photospheric emission and the second the 
infrared excess. If the grains are large (p = 0) and the surface density is determined by PR 
drag (q = 1), then the infrared spectrum, F v oc v~ x . We plot the minimum x 2 fits to the 
IRS photosphere subtracted spectra in blue in Figure 3, assuming that the grains are large 
(p = 0) and list the fitting parameter K 2 in Table 4 along with the minimum reduced x 2 
for the fits. The continuous disk model has a lower reduced x 2 than the single temperature 
black body for one object in our sample (HR 8799); however, this source possess a weak IRS 
excesses (<0.1 Jy) which is detected with a SNR < 5. 

Our simple SED models suggest that PR drag may not be the dominant grain removal 
mechanism in debris disks. To test this hypothesis, we compare the PR drag lifetimes and 
the lifetimes for grains under sublimation, corpulscular stellar wind drag, and collisions for 
all of the sources in our study. The PR drag lifetime of grains in a circular orbit, a distance 
D from a star is 



(Burns et al. 1979). The PR drag lifetimes of average-sized grains (12000 yr - 1.1 Myr; 
see Table 5), estimated using the grain properties in Table 4, are significantly shorter than 
the stellar ages (t age ), suggesting that the grains are replenished through collisions between 
larger bodies. The PR drag lifetime of average-sized grains around HR 8799, the only object 
whose excess spectrum is better modeled by F„ oc A, is not the shortest grain lifetime by a 
factor of 2; therefore, we do not expect that this system should possess a uniform disk. We 
estimate lower limits for the parent body masses around our objects assuming that all of the 
grains are destroyed by PR drag and that the systems are in steady state. If M PB denotes 
the mass in parent bodies, then we may write 



-3+2g+0.5pq-0.5p 



8 




(9) 



M PB > 



4L IR t, 

c 2 



age 



(10) 



-16- 



(Chen & Jura 2001). We estimate the infrared luminosities of the systems from the black 
body fits to the excess. If the grains emit a substantial fraction of their radiation at A >30 
/ill!, then this approximation may not be valid. The inferred parent body masses range 
between 1.3 x 10~ 3 M e for the K1V star HD 53143 and 0.93 M e for the A3IV/V star HR 
1082 (see Table 5). 

Whether PR drag or other processes, such as collisions, are the dominant grain destruc- 
tion mechanism in debris disks is uncertain and depends on the density of dust grains and 
the spatial distribution of dust for each particular object. Numerical models of the dynami- 
cal evolution of dust in collisional equilibrium suggest that Poynting-Robertson drag is the 
primary mechanism for grain transport in disks with L/^/L* < 10~ 6 while radiation pressure 
is the primary mechanism for dust transport in collisionally dominated disks with L IR / 
> 10" 4 (Krivov et al. 2000). Recently, Wyatt (2005) concluded that PR drag is relatively 
unimportant in IRAS -discovered debris disks. He found that the volume density of grains in 
iTMS'-discovered debris disks is so high that grains which migrate inward under Poynting- 
Robertson drag will suffer destructive collisions with other grains and will be rapidly expelled 
by radiation pressure. The destruction of inward-migrating dust grains via mutual collisions 
may explain the presence of central clearings without requiring the presence of icy grains or 
planets in debris disks. 



6.3. Stellar Wind Drag 

For B-type main sequence stars and young solar-like main sequence stars, both of which 
possess strong stellar winds, drag on dust grains produced by loss of angular momentum 
to corpuscular stellar wind may be stronger than that produced by the Poynting-Robertson 
effect. Stellar wind drag may explain the observed ant i- correlation between Spitzer 24 jum 
excess and ROSAT fluxes toward F-type stars in the 3-20 Myr Sco-Cen (Chen et al. 2005) 
and the lack of 12 /iin excesses observed toward nearby, >10 Myr-old, late- type M-dwarfs 
(Plavchan et al. 2005). Recently, Strubbe & Chiang (2006) have reproduced the radial 
brightness profile of the AU Mic disk assuming that collisions between parent bodies on 
circular orbits at 43 AU produce the observed dust grains. Large grains produce a surface 
density, a oc r°, at r < 43 AU, under corpuscular and Poynting-Robertson (CPR) drag 
modified by collisions while small grains that are barely bound under corpuscular stellar 
wind and radiation pressure produce a surface density, a oc r~ 5 / 2 , in the outer disk. 

Stellar winds around B-type stars are produced by the transfer of momentum from 
photons below the photosphere to material outflowing in the wind. The mass loss rates and 
expansion velocities can be measured from radio observations of free-free emission produced 
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in the outer parts of the wind or from Ha emission which may be formed in the inner regions 
of the wind. Monte Carlo simulations of stellar winds that include multiple scattering of 
photons, successfully reproduce stellar wind mass loss rates measured at radio wavelengths 
(Vink et al. 2000). We estimate the stellar mass loss rates (M win d) f° r the B-type stars in 
our sample using a fit to the numerical models of Vink et al. (2000) that depend on the 
stellar luminosity, mass, and effective temperature, assuming that the ratio of the terminal 
wind velocity to the effective wind escape velocity at the stellar surface, fooAw = 1-3. The 
increase in "drag" in the inward drift velocity of dust grains under corpuscular stellar wind 
and Poynting- Robertson drag over that produced by Poynting-Robertson drag alone is given 
approximately by the factor (1 + M wind c 2 / X*) (Jura 2004). 

The mass loss rates due to stellar winds around 14 nearby, solar-like stars have been in- 
ferred via Lyman a absorption, produced when the stellar wind collides with the surrounding 
interstellar medium, producing a hot H I astrosphere with an effective temperature 20,000 - 
40,000 K (Wood et al. 2002, 2005). Wood et al. (2002, 2005) fit the stellar mass loss rate, 
M W i n d, per stellar surface area, A, as a power-law function of x-ray flux per stellar area, 
M W i n d/A oc F X ' 34±0A8 , assuming that the wind speed for solar-like stars is similar to that 
measured for the Sun, v W i n d = 400 km s -1 . The uncertainty in the M win d extrapolation is 
probably a factor of two because the size of the astrosphere and the amount of astrospheric 
absorption scales as the square root of the wind ram pressure; the wind ram pressure, P W i n d 
oc M win dV W ind] and the variation in the solar wind speed is approximately a factor of two. 
The stellar mass loss rate power law dependence on x-ray flux per stellar area saturates at 
F x = 8xl0 5 ergs cm" 2 s" 1 . One possible explanation for the saturation is that stars possess 
more polar spots as they become more magnetically active, indicating changes in the field 
geometry. Changes in the field structure could precipitate changes in the stellar wind (Wood 
et al. 2005). 

We infer M wind from ROSAT fluxes for our sample assuming the Wood et al. (2005) 
power-law, scaling to observations of a Cen (F x = 3.7 x 10 4 erg cm -2 s _1 , M win d/A = 0.9 
Mq/Aq); the astrosphere for this source is well-detected in Lyman a and its M W i n d/A and F x 
lie on the published fit. We list the observed ROSAT fluxes, HR 1 hardness ratios between 
the 0.1 - 0.4 and the 0.5 - 2.0 keV bands, and the angular offsets between the ROSAT catalog 
sources and the FGKM-type stars in our sample in Table 6. The x-ray spectra for stars in our 
sample have flat spectral energy distributions, with a mean hardness ratio -0.4, consistent 
with observations of stars in the solar neighborhood. We estimate X-ray luminosities, Lx, 
using the conversion 1 ROSAT count = (8.31 + 5.30 HR1) x 10~ 12 ergs cm" 2 (Fleming et al. 
1995). The inferred F x > 8xl0 5 ergs cm" 2 s" 1 for 8 solar-like and low-mass stars, suggesting 
that the stellar winds for these stars are probably one or two orders of magnitude smaller 
than inferred from the Wood et al. (2005) relation: r 1 Eri, a For, HD 53143, HD 113766, 
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HD 139664, HD 146897, HD 181327, and HD 191089. We do not make any extrapolations 
for the mass loss rates in these systems. 

For the B-type stars and solar-like and low mass stars with F x < 8xl0 5 ergs cm -2 s -1 , 
M wind c 2 /L* > 1, suggesting that stellar wind drag can not be neglected (see Table 5). The 
lifetime for grains in a circular orbit under the Poynting-Robertson effect and stellar wind 
drag is 

1 

tpR+wind — -y—- : \-tpR (11) 

I M wind c z _|_ ]_ J 

We estimate the combined Poynting-Robertson and stellar wind drag lifetimes using the 
Poynting-Robertson drag time lifetimes and M win dC 2 / values listed in Table 5. For B-type 
stars, M wini ic 2 / L* is typically ~2; therefore, the drag lifetimes of the grains are reduced in 
most cases by ~60%. For example, A Cas and HR 6532 possess M W i n s? / = 1.4 and 1.2, 
respectively. However, for F-type and later stars, M win dC 2 / X* is typically ~100 if the source 
is detected by ROSAT; therefore, the drag lifetimes of the grains are reduced by a factor 
of a couple hundred. For example, HR 506 and HR 6670 may possess M wind c 2 / as high 
as 82 and 190. Unfortunately, ROSAT upper limits do not place stringent upper limits on 
M winc ic 2 1 in cases in which the source is not detected. The mass in parent bodies assuming 
that the system is in steady state and that Poynting-Robertson drag and stellar wind drag 
are the dominant grain removal mechanisms is given by the expression, 

M PB > (^f^ + ij L -ft age (12) 

We estimate the mass in parent bodies using the inferred stellar properties listed in Table 
4. For B-type stars, the stellar wind drag inferred parent body masses are consistent with 
the PR drag inferred parent body masses, M PB <1 M®. However, for solar-like stars, the 
stellar wind drag inferred parent body masses are substantially (20 - 50 x) higher but still 
less than 1 M®. The fact that so little material remains in parent bodies suggests that planet 
formation must be an efficient process if planets have already formed in these systems. 



6.4. Collisions 

If the particle density within the disk is high, then collisions are expected to dominate 
grain destruction by generating small grains that are removed rapidly via radiation and 
corpuscular stellar wind pressure. We estimate the collision lifetime assuming that grains 
are on inclined orbits, such that they encounter the surface density of the disk twice per 
orbit, that tangential collisions are destructive, and that they have a collisional equilibrium 
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size distribution (given in equation 4) 

t coll = 3000 ji{D/^AUyl\M subm jQAM & )-\^— ^/l mm)(M /M,) 1 / 2 (13) 

Backman & Paresce (1993) using grain distances, dust masses (extrapolated for micron-sized 
through centimeter-sized grains), and minimum radii in Table 4. For the majority of our 
systems, the collision lifetimes computed are comparable to or shorter than the Poynting- 
Robertson and/or stellar wind drag lifetimes (see Table 5), suggesting that collisions dom- 
inate the destruction of particles, consistent with published studies (Wyatt 2005; Domink 
& Decin 2003; Najita & Williams 2005). In Figure 4, we plot the sublimation lifetime if 
the grains are icy, the Poynting-Robertson and stellar wind drag lifetime, and the collision 
lifetime as a function of radius around typical B5V, A5V, and F5V stars. In all cases, subli- 
mation is the dominant grain removal process at small distances from the star. For typical 
A5V and F5V stars, the collision lifetime for average-sized grains is shorter than the drag 
lifetime if the disk has a dust mass between 0.001 M® and 1 M®, even if the F5V star has 
a stellar wind with a mass loss rate as high as M wind = 1000 M®. However, for a typical 
B5V star, the Poynting-Robertson and stellar wind drag lifetime may be shorter than the 
collision lifetime, especially at large radii, if the disk has a dust mass, M submm ~ 0.1 M e ; 
however, we infer small M su b mm = l.OxlO -5 M® for some stars in our sample, suggesting 
that Poynting-Robertson and corpuscular stellar wind drag effects should be the dominant 
grain removal mechanism for at least some objects. 

The fractional infrared luminosity of a debris disk for a fixed distance is expected to 
decrease inversely with time, Lm/L* oc l/t age , if collisions are the dominant grain removal 
process and is expected to decrease inversely with time squared, Lj R / oc l/t 2 age , if Poynting- 
Robertson drag is the dominant grain removal process (Domink & Decin 2003). We plot 
Lir/L* as a function of age for the systems in our sample in Figure 6a. For young stars 
which do not appear on the Schaller et al. (1992) isochrones, we assume moving group ages 
rather than our assigned ages. The upper envelope of the distribution can be fitted with 
the function, L IR /L* = (L IR /L*) (t /t age ), where (L IR / L*) t = 0.40 Myr; however, the 
function L IR / = (L IR /L*) (t /t age ) 2 , with (L IR / L^ot 2 = 60 Myr 2 does not produce a 
bad fit. The 1/t 2 trend line includes all of the sources in our sample except r) Crv which has 
a high Li R /L* = 3xl0~ 4 for its age of 1 Gyr. Our data set is consistent with the idea that 
debris disks are collisionally dominated systems but our data set is too small to determine 
whether Poynting-Robertson drag dominates grain removal at any particular age. A l/t age 
time dependence with a decay timescale t Q = 150 Myr has been observed for the MIPS 24 
/im and IRAS 25 /im excess around ~270 A-type stars, consistent with grain destruction via 
collisions (Rieke et al. 2005). Similarly, A l/t age time dependence with a decay timescale t Q 
= 200 Myr has been observed for the 850 fxm excess around 13 nearby, solar-like stars also 
consistent with grain destruction via collisions (Najita & Williams 2005). Our study differs 
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from these in the Rieke et al. (2005) and the Najita & Williams (2005) samples in that (1) 
our sample contains a heterogenous mix of stars with varying spectral type and (2) our data 
set includes the measurement of the excess at more than one wavelength; the measurement 
of the SED allows us to infer the dust temperature and therefore the infrared luminosities 
and dust masses/dust production rates more accurately. 

In a minimum mass solar nebula, 1000 - 3000 km-sized bodies are expected to grow on 
timescales, t P « 15 - 20 Myr (D/30 AU) 3 (Kenyon & Bromley 2004, 2005) and to perturb 
planetesimals in the disk, initiating collisional cascades that produce micron-sized grains. 
This model predicts a t~ a6 to £~ - 35 decay of the fractional infrared luminosity for the whole 
debris disk at ages of 10 Myr to 1 Gyr, depending on the tensile strength of the grains, 
a somewhat more shallow evolution of the fractional infrared luminosity than inferred by 
Dominik & Decin (2003). The Kenyon & Bromley (2004, 2005) models also predict that 
older debris disks systems possess infrared bright rings of dust at larger radii than their 
younger counterparts. Since the dust grains in our sample probably have a > 10 /im, we plot 
the inferred black body dust distance, D, as a function of stellar age to test this hypothesis 
(see Figure 6b). The black body distances, inferred from SED models, are consistent with the 
measured radii of HR 4796 A, Fomalhaut, and Vega from maps of thermal emission from large 
particles at mid-infrared and submillimeter wavelengths to within a factor of 2 (Holland et al. 
2004; Wilner et al. 2002; Jayawardhana et al. 1998). However, the black body distances are 
smaller than the measured radii of (3 Pic and AU Mic in scattered light (Krist et al. 2005; 
Heap et al. 2000) and than the measured radius of Vega, inferred from maps of thermal 
emission of stochasitically heated small grains (Su et al. 2005). We do not find a clear 
correlation between dust grain distance and the age of the central star, in agreement with 
submillimeter studies (Najita & Williams 2005). The addition of submillimeter data for rj 
Crv and r Cet does not appear to improve our fit of grain distance as a function of stellar 
age. In Figure 6b, we overplot grain distances for these objects inferred from submillimeter 
SED models. 



7. Correlations Between Stellar and Dust Properties? 

We searched for correlations between stellar and dust properties. We investigated 
whether the multiplicity, the rotational velocity, or the metallicity of the central star is 
correlated with either the observed fractional infrared luminosity or grain temperature. 

In a binary system, the orbits of dust grains at distances approximately 1.6 to 2.6 times 
the binary separation are expected to be unstable (Artymowicz & Lubow 1994), leading to 
the formation of two populations of circumstellar dust grains. Disks around each component 
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of the binary system are trucated at their outer radii but cool dust may reside at much larger 
distances in circumbinary disks. We searched for a correlation between the single/binary na- 
ture of a stellar system and the inferred circumstellar dust properties (grain temperature and 
fractional infrared luminosity) to determine whether the gravitational effects of a secondary 
star affect disk properties. 

Eleven stars in our sample are binary systems; five of which (45.5%) possess IRS excesses 
detected with good signal-to- noise ratios (SNRs). One of the binary systems with IRS excess 
does not possess spectral features and is better modeled using a single temperature black 
body than a uniform disk whose surface density is given by Poynting-Robertson Drag (HR 
1570). The grain temperature for this systems is T gr = 90.6 K and the fractional infrared 
luminosity for this system is 3.0 xlO -5 . For comparison, forty-nine stars in our sample are 
single systems; twenty-nine of which (59.2%) possess IRS excesses detected with good SNRs. 
Fifteen of the single systems with IRS excesses do not possess spectral features and are 
better modeled using a single temperature black body (HR 333, HR 506, 7 Tri, HR 1082, 
HD 53143, HR 3314, HD 95086, HD 110058, A Boo, HD 139664, HD 146897, HR 6297, HR 
6486, HD 181327, and HD 191089). The mean grain temperature for these systems is T gr = 
103 K with a standard deviation of 19 K. The mean fractional infrared luminosity of these 
systems is 8.4xl0~ 4 , with a standard deviation of 1.5xl0~ 3 . The mean grain temperatures 
and fractional infrared luminosities for single and binary systems are consistent when the 
standard deviation of these quantities are taken into account. 

One difference appears when the single and binary star populations are compared. The 
fraction of infrared excess binary systems which can be modeled by a single temperature black 
body is significantly smaller than the number of infrared excess single star systems. We list 
the angular and physical separation of all the binary systems in our sample in Table 7. One 
possibility for the poor single temperature black body fits associated with binary systems 
is the presence of multiple populations of circumstellar dust grains (e.g., circumprimary, 
circumsecondary, and circumbinary disks). If the binary nature of these systems contributes 
to the complicated structure of the infrared excess then objects which are not well modeled 
by a single temperature black body should possess binary separations 1 - 100 AU. However, 
the separations of one of the two systems (A Cas, HR 6532) which are not well modeled 
by a single temperature black body is too small to truncate a circumstellar disk. HR 6532 
is a spectroscopic binary with a period of 6.8 days. Larger statistics are needed to confirm 
whether the infrared spectra of binary systems is more complicated than a single temperature 
black body. 

We plot the fractional infrared luminosity and grain temperature as a function of binary 
separation (see Figure 7) to examine further the effects of a secondary star. We assume very 
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small separations for Algol eclipsing and spectroscopic binaries (<1 AU). Systems to the 
left of the dotted vertical line have separations <10", small enough that both the primary 
and secondary fall into the LL (and sometimes the SL) slit. Systems to the right of the 
dotted verticle line are too widely separated for the secondary to contribute to either the 
SL or LL spectrum. The average fractional infrared luminosity for systems with separations 
of 10 - 100 AU, 1.6xl0~ 5 , is smaller than that inferred for objects with separations of 100 
- 500 AU, 7.6 xl0~ 4 ; however, this comparison is highly biased by the very high fractional 
infrared luminosity associated with HD 113766, 1.5xl0~ 2 , because our binary sample is so 
small. The average grain temperature for systems with separations of 10 - 100 AU, 120 K, is 
consistent with the peak in grain temperature distribution for the whole sample, but smaller 
than that inferred for systems with separations of 100 - 500 AU, 240 K. Debris disks with 
T gr ~ 250 K are rare (Aumann & Probst 1991). The very high average T gr estimated is 
based on a sample of two unusual objects (HD 113766 and rj Tel). Detailed modeling of the 
10 /im feature observed toward these objects suggests black body grains with T gr > 300 K 
and amorphous olivine with T gr > 200 K (see section 4; Schutz, Meeus, & Sterzik 2005). 

X-ray studies of solar-like stars suggest that magnetic breaking causes stars to rotate 
more slowly as they age, with a v ro t oc 

t -0.6±0.1 t 

ime dependence for stars with ages >0.3 
Gyr (Ayers 1997). We searched for correlations between the inferred circumstellar dust 
properties and the measured projected stellar rotational velocity, vsmi. Using v smi as a 
proxy for stellar age, we expect that both the fractional infrared luminosity and the grain 
temperature will decrease with stellar age or increase with vsmi. We plot the fractional 
infrared luminosity and grain temperature as a function of v sin i in Figure 7, using different 
symbols for B-, A-, F-, and K-type stars. For binary systems, we used the vsmi for the 
primary star. Stars with spectral type earlier than F5V are shown with open symbols while 
stars with spectral type F5V and later are shown with solid symbols. Early-type stars are not 
expected to spin down as they age because they do not possess deep convective envelopes. 
When F5V and later spectral-type stars in our sample are compared with one another, a 
possible decrease in fractional infrared luminosity as a function of v sin i and increase in 
grain temperature as a function of v sin i are seen. However, our sample is predominantly B- 
and A-type stars and only contains 3 mid-F and later type stars with measured dust grain 
properties and measured vsmi] therefore, no strong conclusions about the dependence of 
dust grain properties on v sin i can be drawn in this study. 

Radial velocity studies of main sequence stars that possess giant planets find a corre- 
lation between the presence of an orbiting planet and the metallicity of the central star. 
Spectral synthesis modeling of high resolution visual spectra, sensitive to stellar semiampli- 
tudes > 30 m s _1 and orbital periods shorter than 4 yr, find that fewer than 3% of stars with 
-0.5 < [Fe/H] < 0.0 have Doppler-detected planets while 25% of stars with [Fe/H] > +0.3 
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possess giant planets (Fischer & Valenti 2005). Recent Spitzer MIPS observations suggest 
that nearby planet-bearing, solar-like stars may be more likely to possess 70 /im excesses 
and larger average 70 /zm excesses than stars without known planets (Beichman et al. 2005). 
These excesses are generated by cool dust (T gr < 100 K) located beyond 10 AU, well out- 
side the orbits of the discovered planets. If correlations exist between metallicity and the 
presence of a planet and between the presence of a planet and a 70 /im excess, there might 
also be a correlation between the presence of IRS excess and stellar metallicity; although, no 
correlations between the presence of 70 /im excess and stellar metallicity have been found 
thus far (Bryden et al. 2006; Beichman et al. 2005). 

We plot fractional infrared luminosity and grain temperature versus [Fe/H] for the eleven 
F-type stars in our sample that possess [Fe/H] measurements determined from the Stromgren 
photometry survey of Nordstrom et al. (2004; see Figure 7). For the six objects that do not 
possess IRS excess, we indicate their position along [Fe/H] axis with an upper limit symbol. 
The stars with and without detected IRS excess have similar metallicities, consistent with 
a stochastic origin for the small infrared emitting grains in debris disks. We do not include 
B- and A-type stars in this plot because they possess shallow surface convective zones. The 
observed metallicity of these objects may be more easily distorted by pollution than their 
F-type counterparts. We observe a large dispersion in both the fractional infrared luminosity 
and grain temperature even though our sample only contains 5 F-type stars with measured 
grain properties and measured [Fe/H]. More objects from our sample may be folded in when 
grain properties inferred from longer wavelength observations are included; for example, r 
Cet, with [Fe/H] = -0.47, log{L IR /L*) = -4.6 (Decin et al. 2003) and T gr = 60 K (Greaves 
et al. 2004) inferred from far-infrared and submillimeter photometry. 



8. Gas Mass Upper Limits 

Atomic and molecular gases may affect the dynamics of circumstellar dust grains. Nu- 
merical models of disks with gas:dust ratios of 0.1 - 10 suggest that gas-grain interactions 
may generate the observed infrared bright rings by concentrating small grains, with radii just 
above the blowout radius, at the outer edge of the gas disk (Takeuchi & Artymowicz 2001). 
In addition, measurement of gas masses in a large sample of disks may help determine the 
gas dissipation timescale and help constrain models for giant planet formation (Hollenbach 
et al. 2005). We searched for H 2 S(0), S(l) and S I emission at 28.2 /im, 17.0 /im, and 
25.2 fim, respectively, to constrain the bulk gas mass in debris disks. We do not detect 
any of these emission lines in any of our SH and LH module data. We estimate 3a upper 
limits to the H 2 S(0), S(l) and S I line fluxes from our nod difference spectra (see Table 8). 
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Since the amplitude of each difference spectrum varied from pixel to pixel near each line, we 
averaged each difference spectrum in a region ±0.5 /im around each line to determine the 
uncertainty in the line flux. We estimate H 2 mass upper limits from the S(0) line flux upper 
limits, assuming that the source is unresolved and that the gas has a temperature, T ex = 
50 K and fOO K, expected from bulk gas that is co-spatial with the infrared-emitting dust. 
Since the most constraining gas mass upper limits are typically <100 M e and the measured 
submillimeter dust masses are typically <1 M e , we can not constrain the gas:dust ratio in 
debris disks well. 

Molecular gas has been detected toward the nearby (70 pc away from the Sun) ~20 Myr 
old star 49 Cet. Submillimeter observations find CO J = 3^2 and J = 2— >1 emission at the 
radial velocity of the star (Dent et al. 2005; Zuckerman et al. f995). Models of the 49 Cet 
doubly peaked CO J = 3^2 profile are consistent with a compact disk with an outer radius, 
D out = 17 AU, and a disk inclination, i = 16°; the high velocity wings of the profile (~4 km 
sec -1 relative to the star) suggest that gas is present at <5 AU from the central star (Dent 
et al. 2005). ISO observations detect H2 S(0) and S(l) emission with a line flux 6.6±2.0 x 
f 0~ 14 ergs sec -1 cm -2 and place 3<r upper limits on the S(l) line flux, <0.8 x f 0~ 14 ergs sec -1 
cm~ 2 , corresponding to a gas:dust ratio >f00 (Thi et al. 200f). Our 49 Cet 3cr H 2 S(0) line 
flux upper limit apparently conflicts with the Thi et al. (200f ) results because our S(0) upper 
limit is a factor ~9 lower than the reported detection. Our upper limit is consistent with the 
mass expected in H 2 inferred from 800 /im - If 00 /im continuum measurements of the dust 
mass, M suhmm = 0.02 - f M e , if the disk has an interstellar gas:dust ratio (Zuckerman et al. 
f 995) and is also consistent with the inferred CO masses from submillimeter measurements. 
Zuckerman et al. (f995) and Dent et al. (2005) estimate that between 2 and 30 M ffi H2 exists 
in the disk, assuming a gas excitation temperature, T ex = 50 - 60 K, and an interstellar 
H 2 :CO number ratio between f - 2xf0 4 , significantly smaller than the <6000 M e and <80 
M e that we infer assuming T ex = 50 K and fOO K, respectively. 

The non-detection of H2 emission at mid-infrared wavelengths using Spitzer IRS may be 
somewhat constraining for at least one object in our sample. The expected H2 emission line 
luminosity may be estimated from our inferred parent body masses (using the estimates in 
Table 5) assuming an interstellar gas:dust ratio and a gas excitation temperature, T ex = T gr . 
If the system is a point source, then the total luminosity produced by N(H 2 ) H 2 molecules 
is 

huN(H 2 ) Xu A ul 

F = W (14) 

where E — hv is the energy of the radiated photons, Xu is the fraction of H2 in level u, 
and A u \ is the transition probability. HD f 13766 is a F3/F5V member of Lower Centaurus 
Crux in Scorpius-Centaurus, with an estimated age of 16 Myr and a noteably large L/^/L* 
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= 0.015 (Chen et al. 2005). In our high-resolution mode sample, HD 113766 possesses the 
largest parent body mass (0.1 M e or 260 times that mass in the main asteroid belt) and 
the hottest single black body grain temperature, T gr = 330 K, making it the most likely 
object to possess H 2 S(0) and S(l) emission. The approximation T ex = T gr is valid for 
optically thick, flared disks for Ay > 0.1 (Kamp & Dullemond 2004) but invalid for debris 
disks around A-type stars where the disk gas and dust are too tenuous to heat the gas via 
gas-grain collisions (Kamp & van Zadelhoff 2004). For example, modeling of the gas around 
HR 4796A, an 8 Myr old A-type member of the TW Hydrae Association, located 67 pc away 
from the Sun, suggests that the gas, T ex = 65 K, is substantially cooler than the dust, T gr 
= 110 K (Chen & Kamp 2004). With these caveats, we estimate that the HD 113766 disk 
should possess H 2 S(0) and S(l) line fluxes of 2.8xl0~ 15 ergs s" 1 cm" 2 and 2.8xl0" 13 ergs 
s _1 cm~ 2 , respectively, corresponding to an S(l) emission line flux ~3x higher than our 
observed 3a upper limit. 

Gorti & Hollenbach (2004) have argued that [S I] emission at 25.2 /im may be an ex- 
cellent tracer for bulk gas in proto-planetary disks. Their model spectra of disks around 
intermediate-aged (~10 7 years) stars predicts strong [S I] emission, with line: continuum ra- 
tios >5% in Spitzer IRS high- resolution mode (R ~ 600), assuming that the disk sulfur 
abundance is similar to that observed in the diffuse interstellar medium (2.8xl0~ 5 ) because 
co-spatial dust grains are too warm for volatiles to condense on them. However, the gas- 
phase sulfur abundance in intermediate-aged disks may be depleted relative to the diffuse 
interstellar medium. The previously-identified 23.5 /im wustite (FeO) emission feature de- 
tected in ISO spectra of Herbig Ae stars may not be wustite but pyrrhotite ([Fe,Ni]!_ x S). 
Keller et al. (2002) compare the intrinsic strength of the 18 /im silicate feature to that of the 
23.5 /im pyrrhotite feature in the ISO spectrum of HD 163496 and measure a silicon:sulfur 
ratio of 0.63, consistent with a solar abundance (0.52). They conclude that most if not all 
of the sulfur around HD 163296 resides in solid FeS grains. If the bulk of the circumstellar 
sulfur is located in grains, then [S I] emission at 25.2 /im (or lack thereof) may not be a good 
bulk gas tracer. 



9. Conclusions 

We have obtained Spitzer Space Telescope IRS spectra of 59 main sequence stars with 
spectral types B- through M-type and ages between 1 Myr and a few Gyr, that possess IRAS 
60 /im excess. We find that: 

1. The majority of observed debris disks do not possess spectral features, suggesting 
that the grains are too cold and/or too large (a > 10 /im) to produce spectral features. 
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Detailed modeling of objects with spectral features requires the presence of large, warm 
amorphous silicates with T gr = 290 - 600 K, in addition to cool black body grains with T gr 
= 80 - 200 K, and the presence of crystalline silicate mass ratios 0% - 76%. 

2. The IRS spectra of debris disks (without spectral features) are generally better fit 
using a single temperature black body than with a uniform disk. Stellar radiation pressure 
(in collisionally dominated systems), sublimation if the grains are icy, gas drag, and/or the 
presence of a perturbing body may contribute to the presence of inner holes in these disks. 

3. The peak in the distribution of estimated black body grain temperatures, T gr = 110 - 
120 K, suggesting that sublimation of icy grains may produce the central clearings observed. 

4. Since the parent body masses typically are less than the mass of the Earth, it appears 
that planet formation efficiently consumes most of the mass of the primordial disk. 

We would like to thank K. Uchida for his assistance with data reduction and our anony- 
mous referee, S. Kenyon, E. Mamajek, J. Mould, I. Song, A. Speck, K. Stapelfeldt, and S. 
Strom for their helpful comments and suggestions. This material is based upon work sup- 
ported by the National Aeronautics and Space Administration under Award No. NAS7-1407 
and the California Institute of Technology. 

A. Selected Source Notes 

HR 506 may possess a planet with a period of 1040 days, a minimum (M sin i) mass of 
0.91 M Jup , and a doppler velocity semi-amplitude of 18 m s _1 ; the inferred eccentricity and 
semi-major axis for the planetary orbit are e = 0.18 and a = 2.1 AU (Mayor et al. 2003). 
If HR 506 possesses black body grains at a distance D = 2.1 AU from the star, then they 
are expected to have a grain temperature T gr = 220 K, assuming a stellar radius R* = 1.2 
R Q and an effective temperature T e ff = 6140 K, consistent with a stellar spectral type of 
F8V. However, the IRS excess continuum emission observed toward this star is better fit by 
a black body with a cooler grain temperature T gr = 70 K, suggesting that planetesimals at 
similar distances have already been removed. The black body distance for T gr = 70 K grains 
is D = 21 AU. Therefore, the observed infrared emitting dust and parent body population 
are located at significantly larger distances than the radial velocity planet. A second planet 
in this system is probably needed to dynamically stir the parent bodies because the observed 
radial velocity planet is located too far away to generate strong orbital resonances (Bryden 
et al. 2000; Beichman et al. 2006). 

t Cet is a 7.2 Gyr old G8V star (Lachaume et al. 1999) that possesses IRAS and JCMT 
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SCUBA excesses that are well fit by a T gr = 60 K modified black body (Greaves et al. 2004). 
However, the thermal emission from this cold dust population is not well detected at IRS 
wavelengths (A < 35 /im). The slope of the 20 /im - 35 /im IRS spectrum is expected to be 
~4% higher assuming the Greaves et al. (2004) model. The difficulty in detecting thermal 
emission from this population of circumstellar grains suggests that the IRS is sensitive to 
grains with dust temperature, T gr > 60 K. 

HR 3927 is a A0V star with an estimated age of 50 Myr. Its Spitzer IRS spectrum 
possesses a 10 /im feature that may be modeled using solid amorphous olivine and crystalline 
forsterite grains with a = 3.1 and 8 /im, respectively, larger than the blow-out size for both 
species, T gr = 290 K, and a crystalline silicate fraction of 38%. The residual continuum 
emission may be modeled using black body grains with a lower grain temperature, T gr = 
80 K. The presence of a cooler black body continuum in addition to the hot silicate grains 
suggests that this system possess two planetesimal belts in analogy with the asteroid and 
Kuiper belts in our solar system. The IRAS (not color-corrected) 60 /im flux of 0.69 Jy 
is somewhat higher than the 0.34 Jy inferred from our model; however, the presence of a 
second source in the beam can not be ruled out without Spitzer MIPS images. If HR 3927 
does possess a larger 60 /im excess than predicted by our model, then it may also possess an 
additional (or third) population of cold dust with T gr < 60 K. 

1] Crv is a F2V star with an estimate age of 1 Gyr based on its x-ray activity. This 
system may possess three populations of dusty debris. The coldest population with T gr = 
40 ±5 K has been detected at sub millimeter wavelengths where the system is resolved at all 
position angles at 450 /im with an elongation at a position angle of 130°±10° and at 850 
/im with a radius of 100 AU (Wyatt et al. 2005). A second component is detected in IRAS 
photometry; self-consistent modeling of the IRAS and sub millimeter photometry suggest 
that this population possesses a grain temperature T gr = 370±60 K (Wyatt et al. 2005). 
Detailed fits to the Spitzer IRS 5.5 - 35 /im spectra suggest that this emission is produced 
by amorphous olivine, and crystalline forsterite and enstatite features with T gr = 360 K and 
a crystalline silicate fraction of 31%, and black body grains with a temperature T gr = 120 K. 
The detection of three planetesimal belts may require the presence of at least two perturbing 
planets because three belts are unlikely to be located at the resonances of a single planet. 
Detailed dynamical models of this system are required to determine whether rj Crv possesses 
multiple planets. 

HD 113766 is a binary member of Lower Centaurus Crux in the Sco-Cen OB Associ- 
ation with an estimated age of ~16 Myr (Mamajek et al. 2002). The secondary lies 1.3" (or 
170 AU) away from the F3/F5 primary. HD 113766 possesses an extremely high fractional 
infrared luminosity L/^/L* = 1.5xl0~ 2 and MIPS 24 /im and 70 /im fluxes that are well-fit 
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with a T gr = 330 K black body, suggestive of debris at terrestrial planet temperatures (Chen 
et al. 2005). Published ground-based ESO TIMMI2 spectra suggest that the 10 /im feature 
is dominated by crystalline silicate (forsterite) and large, amorphous silicates; SiC>2, which 
is correlated with the presence of forsterite in Herbig AeBe ISO spectra, is not detected 
(Schutz et al. 2005). We use crystalline forsterite in addition to amorphous carbon, olivine, 
and pyroxene, and a single temperature black body to fit not just the 10 /im feature but 
also the 5.5 /im - 35 /im IRS spectrum and infer a crystalline silicate fraction of 4.1%. The 
presence of a cooler black body continuum (T gr = 200 K) in addition to the hot silicate 
and carbon grains (T gr = 600 K) suggests that this system possess two planetesimal belts in 
analogy with the asteroid and Kuiper belts in our solar system. The lack of spectral features 
due to sub-micron-sized grains can be explained if small particles are radiatively driven from 
the system by radiation pressure. 

ji Ara possess a planet with a period of 743 days, a minimum (Msmi) mass of 1.97 
Mj up , and a doppler velocity semi- amplitude of 54 m s -1 ; the inferred eccentricity and semi- 
major axis for the planetary orbit are e = 0.62 and a = 1.65 AU (Butler et al. 2001). If 
the system possessed black body grains at a distance D = 1.65 AU from the star, then 
they would have a grain temperature T gr = 220 K, assuming a stellar radius = 0.99 
R Q and an effective temperature T e ff = 5830 K, consistent with a stellar spectral type of 
G3V. Grains with such high temperatures are expected to produce a black body continuum 
that peaks near 25 /im in the center of the IRS bandpass, similar to that observed toward ( 
Lep. However, no excess emission has been detected toward \i Ara at IRS wavelengths. One 
possible explanation for the lack of infrared excess is that this system has already destroyed 
any planetesimal belt it may have possessed by its estimated age of 6.2 Gyr (Lachaume et 
al. 1999). 

HR 7012 is a A5IV-V member of the (3 Pic moving group with an estimated age of 12 
Myr (Zuckerman et al. 2001). Its Spitzer IRS spectrum possesses a strong 10 /im silicate 
feature that may be modeled using large amorphous olivine and pyroxene grains with a = 
5.0 and 1.1 /im, respectively, and T gr = 520 K and sub-micron-sized crystalline enstatite and 
cristobalite. Cristobalite is used to fit emission in the shoulder of the silicate feature at A < 
8.7 /im, however other materials such as opal (hydrated silica) may also provide acceptable 
fits. We estimate that this system possess a crystalline silicate fraction of 76% however our 
model produces extra emission at 10.5 and 19.5 /im that may indicate that the fraction of 
crystalline pyroxene used is too high. The presence of small sub-/im-sized grains inferred 
from our fit to the HR 7012 spectrum suggest that this system may have experienced a recent 
collision. 

7] Tel is a binary member of the (3 Pic moving group with an estimated age of 12 Myr 
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(Zuckerman et al. 2001). The M7/8V secondary located at a distance of 4" (or 200 AU) 
away from the AOV primary was detected via common proper motion studies (Lowrance et 
al. 2000). The 77 Tel Spitzer IRS spectrum possesses a 10 /im feature that may be modeled 
using solid amorphous olivine grains with a = 3.0 /zm, larger than the blow-out size, and T gr 
= 370 K. The residual continuum emission may be modeled using black body grains with a 
lower temperature, T gr = 115 K. The presence of a cooler black body continuum in addition 
to the hot silicate grains suggests that this system possess two planetesimal belts in analogy 
with the asteroid and Kuiper belts in our solar system. This spectrum does not appear to 
possess crystalline forsterite or enstatite features despite the young age of the system. For 
comparison, HR 7012 (another A- type member of the (3 Pic moving group) has a crystalline 
silicate fraction of 76%. 
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Table I. Stellar Properties 
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Table 2. IRS Photometry in 8.5-13 and 30-34 /zm Bands 
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iyy 


a n 
-4.y 


n nop 
-U.UZO 


Q t; K/l 01 p 
O004ol0 


r hiTi 


1Z fU±l 


i o ^n 
IZoU 


on q 

zy.o 


n no/i 
U.UZ4 


i np i i 

iyo±i 


1 QO 

loZ 


CO 

00.0 


n /i 77 
U.4/ / 


/i noi a n 
4yolo4U 


a ror + 


o /) 1 n_i_ /i 
o41U±4 


o4o(J 


on >-7 
-OO. ( 


n m 1 
-U.U11 


/I 1 0_1_£ 

41o±5 


to 


OO 
ZO.O 


n n7/i 
U.U/4 


orr r n7n 
OOOOU / Z 


HR 1082 

XXXL 1 uUZj 


1 254-D 3 


123 


1 7 


01 4 


1 1 7+2 


1 2 9 


1 04 

X Wt: 


7 96 


3555328 


7 Dor* 


1540±2 


1580 


-36.2 


-0.023 


194±2 


170 


24.8 


0.147 


3555584 


HR 1570 


635±0.9 


663 


-29.4 


-0.044 


175±2 


69.0 


105 


1.50 


3555840 


HR 1686* 


889±1 


936 


-46.9 


-0.050 


106 ±2 


101 


6.1 


0.061 


3556096 


k Lep* 


521±0.6 


517 


4.1 


0.008 


73.1±1.3 


53.1 


19.2 


0.356 


3577600 


32 Ori* 


468±0.5 


492 


-23.6 


-0.048 


50.4±1.7 


47.1 


-0.6 


-0.012 


3556352 


C Lep 


1950±0.6 


1710 


232 


0.136 


886±4 


182 


704 


3.87 


4932864 


5 Dor* 


1210±1 


1240 


-33.1 


-0.027 


135±2 


131 


3.2 


0.024 


3556608 



Table 2 — Continued 







8.5 - 13 


/im 






30 - 34 










Measured 


Predicted 






Measured 


Predicted 








Name 


Flux 


Flux 


Excess 


fractional 


Flux 


Flux 


Excess 


fractional 


bpitzer 




(mJy) 


(mJy) 


(mJy) 


Excess 


(mJy) 


(mJy) 


(mJy) 


Excess 


AOR Kej 


HR 2124* 

111 1, JL — _L 


1260+2 


1 350 

1 • ) (J VJ 


-95 1 

<J?J . 1 


-0 070 


203+5 


144 


54.2 


365 

VJ • ' ) VJ '. ' 


3556864 


XZ Fen* 


48 1 +0 3 

tO . J — l_VJ . *J 


45 3 

1 *. J . *J 


2 7 


061 

W . WVJ 1 


5 5+2 1 


4 6 

t: . VJ 


8 

VJ . (J 


1 68 

VJ . 1 VJO 


3585536 


Aur* 


1010+1 


1045 


-35 8 

• ) / . VJ 


-0 034 

VJ • VJ - Jt: 


140+2 


111 


28.4 


255 


35571 20 


Illy ) •) 1. T: •_) 


380+0 4 


395 

"J • J '..J 


-14 6 

It: . VJ 


-0 037 


82.0+1.1 


41 6 

T: 1 . VJ 


39 3 

" J • 7 . • ) 


919 


3557632 


HR 3220* 

111L <J^AiU 


1260+2 


1260 


-7 1 

1 . 1 


-0 006 

VJ • \)\)\) 


143+2 


137 

1 "J I 


7 7 


057 

VJ . VJtJ t 


3557888 


HR 3314 

111L '.)'.) _L T: 


qc7_|_n c 


978 

>J 1 VJ 


9 5 

U • VJ 


010 

VJ . VJ 1 vj 


281+2 


103 

1 VJ "J 


177 

1 1 1 


1 71 


3558144 

■J '.J '.J <-J I'll 


^ Vel* 

C V CI 


VJO*JVJ_l_ 1 


681 

VJ(J 1 w 


35 1 


005 


766+4 

I VJVJ _I_t: 


721 

1 — 1 


44 6 

tit: . VJ 


062 


493081 6 

t:<JOVJ(J 1 VJ 


FI Vir* 


343+0 4 

(J TT-.J _l — V . T: 


348 


-5 6 

tj . V 


-0 016 

VJ . VJ 1VJ 


37 5+1 4 

•J 1 . ' J _l — L . T: 


41 5 

T: 1 . tj 


-2.2 


-0 054 

VJ . V ) 'J T: 


3559168 


HR 3862* 

1111 UUUZj 


1310+1 


1370 

1 "J 1 V 


-66 

VJ VJ . VJ 


-0 048 

vj .u^u 


140+2 


146 

1 TV.' 


-7 1 

1 • 1 


-0 048 

VJ . UtU 


3558400 


HR 3927 

111 L kjU t-i \ 


1 95+f) 4 

1 iJ 5j -J— \J . T: 


177 

1 1 1 


18 1 

1 VJ . 1 


102 


187+2 

1VJ I 


18 8 


168 

luU 


8 94 

VJ . i_/t: 


3558656 


HD 95086 

1 1 1^ i7tJWO\J 


7c; Q_i_n q 


69 

VJ <J . VJ 


1 2 5 


199 


1 1 4+2 

1 It: _LZj 


7 6 

1 . VJ 


107 

lu 1 


1 4 6 

It:. VJ 


355891 2 


/3 UMa 


4110+5 


4160 


-48.8 


-0.012 


801+3 


440 


361 


0.820 


4930304 


(3 Leo 


6390+7 


6370 


18.3 


0.003 


1410+6 


676 


736 


1.09 


4929793 


G Cent 


259+0.5 


268 


-7.9 


-0.029 


34.7+2.2 


27.3 


7.1 


0.269 


3579392 


i] Crv 


1890+3 


1610 


276 


0.171 


476+1 


179 


197 


1.10 


3559424 


HD 110058 


55.3+0.3 


40.3 


15.0 


0.372 


380+2 


4.3 


376 


88.1 


3579648 


HD 113766 


1870+2 


86.2 


1780 


20.650 


988+4 


9.6 


978 


102 


3579904 


HR 5236* 


931+0.5 


901 


28.9 


0.032 


105+0.9 


96.1 


7.6 


0.078 


3559680 


A Boo 


936+0.9 


937 


-0.5 


-0.001 


272+2 


100 


173 


1.75 


3559936 



Table 2 — Continued 







8.5 - 13 


/zm 






on a 

30 - 34 


/im 








Measured 


Predicted 






Measured 


T N 1 * J 1 

Predicted 








Name 


Til 

Flux 


Til 

Flux 


Excess 


fractional 


Flux 


Flux 


Excess 


fractional 


bpitzer 




(mJy) 


(mJy) 


(mJy) 


Excess 


(mJy) 


(mJy) 


(mJy) 


Excess 


AOR Kej 




1510±1 


1540 

1 "J t: V 


-29 8 

Zj<X . \J 


-0 019 


172±2 


164 


7 4 


045 


3560192 

UtJUUl iXZj 


<v Tra* 

j 11a 


Qnqn_i_9 

%J\J U \J 1 Zj 


3190 


-99 8 


-0 031 


KJ'O tX_l_*J 


345 


9 2 

1/ . Zj 


026 


3^60448 


a CrB 


4860+2 


4830 


30 

■XV . \J 


006 


988+4 


51 1 

ull 


476 


932 

\J . tXOZj 


4929280 


HD 13Q664 

illy ItJ JuUt: 


1260±1 


1250 

1 ZjtX 


7 7 


006 


223+2 


135 


89 2 

W tX . Zj 


666 

W . V V vx 


3560704 


41 Ser* 


3160+3 

"X X \X VJ 


3320 


-159 


-0 048 


364+4 


368 


-5 4 

tX . t: 


-0 015 

W.U ItJ 


3561216 

"X 'X V X Zj X \X 


HD 146897 


34.7±0.3 


27.0 


7.6 


0.283 


329+3 


2.9 


327 


113 


3581696 

U vU X \X iX\X 


/T H PT? 
(X 


1 5j W—l—W . (X 


797 


-46 9 


-0 059 


87 0+3 5 


85 2 

(X*J . Zj 


2 9 

Zj . U 


034 

W . V "X 1 


3561472 

OtXVJ It: 1 Zj 


HR 6297 

1111/ \J i-i u \ 


_ 1 J — i—\J . • ) 


287 

Zjtx 1 


-16 8 

X \X . \J 


-0 058 


80 4+2 2 


31 1 

O X . X 


49 9 

t: tX . tX 


1 64 


3561728 

"X '.J \J X 1 ZjO 


HR 6486 

111 I UtiUU 


1470±2 


1520 


-59 9 

(X<X . <X 


-0 039 


310+5 

*X lU-LtJ 


168 

X uu 


142 


846 


3562240 

O ?X \X Zj i-l±AJ 


HR 6532 

111L V .' ' J "X J-' 


i ofi+0 3 

1U\J _LU . "-7 


102 


4 


039 


40 2+2 

t: VJ . Zj _l Zj . W 


10 5 

X \J . tX 


29.4 


2 74 

Zj . f t: 


3582720 

UUUZj 1 Zj\X 


78 Her 

I (X 11CI 


onfi-Uf) 4 

Zj \J O _I_»X . j: 


208 


9 

W . .x 


004 


47 9+1 2 7 

T 1 . J — 1 — L Zj . 1 


21 9 

Zj X . <_/ 


25 9 

ZjtJ . 1/ 


1 18 


3562496 


u Ara^ 


1400±1 


1490 


-91.2 


-0.061 


146+3 


172 


-14.9 


-0.093 


3562753 


7 Oph 


1340±2 


1320 


19.1 


0.014 


729+3 


134 


589 


4.23 


4931328 


HR 6670 


442±0.5 


456 


-13.5 


-0.030 


78.2+1.9 


47.7 


29.4 


0.602 


3563008 


HR 7012 


1410±1 


658 


753 


1.144 


723+5 


73.3 


649 


8.85 


3563264 


r\ Tel 


468±0.7 


333 


135 


0.405 


469+2 


35.1 


433 


12.3 


3563776 


HD 181327 


161±0.3 


154 


7.1 


0.046 


555+2 


15.8 


539 


32.6 


3564032 


HD 191089 


134±0.4 


131 


3.0 


0.023 


363+4 


14.2 


349 


24.9 


3564288 


HD 200800 


47.2±0.3 


43.0 


4.2 


0.097 


9.5+5.9 


4.6 


5.0 


1.09 


3582976 



Table 2 — Continued 







o.O - io 








OU - 04t 


/iin 








Measured 


Predicted 






Measured 


Predicted 








Name 


Flux 


Flux 


Excess 


Fractional 


Flux 


Flux 


Excess 


Fractional 


Spitzer 




(mJy) 


(mJy) 


(mJy) 


Excess 


(mJy) 


(mJy) 


(mJy) 


Excess 


AOR Key 


HR 8799 


290±0.4 


280 


10.5 


0.038 


72.1±2.1 


30.1 


42.3 


1.42 


3565568 


HD 221354* 


441±0.5 


469 


-28.2 


-0.060 


43.4±1.1 


49.4 


-7.1 


-0.140 


3565824 



*stars without strong IRS excess 
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Table 3. Spectral Feature Fitting Parameters 



HR 3927 rj Crv HD 113766 HR 7012 7/ Tel 



Black Body Components 



T m (K) 


80 


120 200 


200 


115 


Q bbl (10~ 16 sr) 


350 


80 68 


42 


160 


T bb2 (K) 








370 


Q bbl (10- 16 sr) 








1.1 


Temperatures 


T silicate (K) 


290 


360 600 


520 


370 


T carbon (K) 




600 


520 




Silicate a mi „ i0 (/im) 


3.1 


0.54 1.4 


1.1 


2.4 


Carbon a m j„ i0 (/im) 




1.9 


1.4 




Silica a mirij0 (/im) 






1.6 




Amorphous Olivine Properties 


Composition 


MgFeSi0 4 


MgFeSi0 4 MgFeSi0 4 


MgFeSi0 4 


MgFeSi0 4 


Shape 


spheres 


spheres spheres 


spheres 


spheres 


a (/im) 


3.1 


3.5 1.5 


5.0 


3.0 


/ (V vac /V tot ) 





0.35 


0.6 





m (10 20 g) 


24 


7.8 82 


4.4 


8.0 


Amorphous Pyroxene Properties 


Composition 




Mg . 5 Fe . 5 Si0 3 


Mg . 8 Feo.2Si03 




Shape 




spheres 


spheres 




a (/xm) 




1.5 


1.1 




/ (V vac /V tot ) 












m (10 20 g) 




97 


2.4 
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Table 3 — Continued 




HR 3927 77 Crv HD 113766 


HR 7012 r] Tel 


Forsterite Properties 


Composition 


M gl . 9 Feo.iSi0 4 Mg L9 Feo.iSi0 4 Mg 2 Si0 4 




Shape 


spheres spheres nonspheres 




a (jum) 


8 8 sub-/xm 




/ (V vac /V tot ) 







m (10 20 g) 


15 0.47 7.7 




Enstatite Properties 


Composition 


MgSi0 3 


Mg . 7 Feo.3Si04 ... 


Shape 


spheres 


nonspheres 


a (/im) 


1 


sub-/xm 


f (V vac /V tot ) 


0.4 





m (10 20 g) 


3 


20 


Silica Properties 


Composition 




Cristobalite 


Shape 




CDE2 


a (//m) 




Rayleigh-limit 


/ (V vac /V tot ) 







m (10 20 g) 




1.5 


Amorphous Carbon Properties 


a (/im) 


1.9 


1.5 


/ (V vac /V tot ) 








m (10 20 g) 


55 


0.7 
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Table 4. Dust Properties Inferred from Excess Continua 



Black Body Uniform Disk 

Name M* T gr Lm/L* D < a > Mdust MiQ Cm Minimum K2 Minimum 

(L ) (M ) (K) (AU) (H (M e ) (M e ) X 2 (Jy/^m) X 2 



A Cas 


250 


4.0 


120±100 


4.4xl0~ 6 


72 


25 


8.9xl0~ 5 


0.007 


4.4 


4.4 


5.4 


HR 333 


41 


2.5 


110±15 


9.8xl0~ 5 


34 


6.4 


l.lxlO" 4 


0.02 


7.8 


33 


31 


49 Cet 


19 


2.2 


118±6 


3.0xl0~ 4 


22 


3.5 


7.7xl0~ 5 


0.01 


7.6 


100 


19 


HR 506 


1.7 


1.2 


70±110 


1.5xl0" 4 


21 


0.6 


5.1xl0~ 6 


0.006 


9.9 


18 


34 


7 Tri 


54 


2.6 


120±15 


2.2xl0~ 5 


31 


8.1 


2.7xl0~ 5 


0.003 


8.4 


56 


24 


t 3 Eri 


13 


2.0 


150±20 


2.1xl0~ 5 


12 


2.5 


1.2xl0 -6 


0.0002 


15 


24 


17 


HR 1082 


11 


1.9 


70±70 


4.2xl0~ 4 


52 


2.1 


3.8X10" 4 


0.07 


2.9 


20 


48 


HR 1570 


22 


2.2 


90±60 


2.9xl0~ 5 


42 


3.9 


3.2xl0- 5 


0.005 


25 


15 


46 


C Lep 


30 


2.3 


191±3 


6.5xl0~ 5 


10 


5.0 


5.6xl0~ 6 


0.0008 


36 


250 


83 


HD 53143 


0.48 


0.8 


120±60 


l.OxlO -4 


4 


0.2 


6.1xl0- 8 


0.00003 


2.4 


6.4 


23 


HR 3314 


43 


2.5 


130±15 


2.5xl0~ 5 


27 


6.8 


1.9xl0 -5 


0.002 


12 


46.5 


33 


HD 95086 


7.5 


1.7 


80±30 


6.4xl0- 4 


30 


1.8 


1.7xl0 -4 


0.03 


5.1 


23 


46 


P UMa 


80 


2.7 


110±30 


l.lxlO -5 


53 


12 


5.3xl0- 5 


0.005 


30 


72 


72 


(3 Leo 


14 


2.0 


120±15 


2.7xl0^ 5 


19 


2.7 


4.2xl0~ 6 


0.0007 


2.3 


190 


38 


HD 110058 


10 


? 


112±7 


l.OxlO- 3 


20 


? 


? 


? 


4.3 


57 


170 


A Boo 


21 


2.2 


120±20 


3.2xl0- 5 


23 


3.8 


9.7xl0- 6 


0.001 


11 


45 


23 


a CrB 


85 


2.7 


139±7 


1.3xl0 -5 


33 


12 


2.7X10" 5 


0.002 


3.5 


140 


54 


HD 139664 


3.6 


2.7 


100±60 


3.8xl0" 5 


15 


1.0 


1.2xl0 -6 


0.0003 


7.0 


23 


7.5 


HD 146897 


4.6 


1.5 


100±8 


5.4xl0- 3 


17 


1.2 


2.9xl0~ 4 


0.06 


4.2 


70 


240 


HR 6297 


11 


1.8 


110±40 


5.1xl0- 5 


21 


2.3 


8.1xl0- 6 


0.001 


8.4 


10 


15 


HR 6486 


6.8 


1.7 


120±20 


6.7xl0- 5 


14 


1.6 


3.3xl0- 6 


0.0006 


4 


38 


6.9 


HR 6532 


93 


3.0 


120±60 


2.9xl0- 5 


47 


12 


1.2xl0- 4 


0.01 


3.1 


7.2 


3.6 



Table 4 — Continued 



Black Body Uniform Disk 

Name M* T gr Lm/L* D < a > M^ust M\Q Cm Minimum K2 Minimum 

(L ) (M ) (K) (AU) H (M e ) (M e ) X 2 (Jy/^m) X 2 



78 Her 


51 


2.6 


140±160 


1.6x10" 


-5 


23 


19 


2.5x10- 


-5 


0.002 


0.5 


6.5 


0.7 


7 Oph 


39 


2.5 


124±9 


5.2x10" 


-5 


27 


6.1 


3.8x10- 


-5 


0.005 


29 


130 


120 


HR 6670 


4.1 


1.5 


110±120 


2.4x10- 


-5 


13 


1.1 


7.4x10- 


-7 


0.0002 


7.3 


2.4 


9.6 


HD 181327 


2.5 


1.9 


81±7 


3.1x10- 


-3 


20 


0.9 


1.7x10- 


-4 


0.04 


21 


104 


1500 


HD 191089 


3.1 


1.8 


99±9 


1.5x10- 


-3 


14 


0.9 


4.4x10- 


-5 


0.01 


3.6 


75 


220 


HR 8799 


6.1 


1.6 


150±30 


4.9x10- 


-5 


8 


1.5 


7.2x10- 


-7 


0.0002 


4.7 


13 


3.1 



Table 5. Dust Removal Mechanisms 



PR Drag CPR Drag Collision Sublimation 

Name tpR M PB t PR+wind M PB t coU t suU 

(Myr) (M e ) (Myr) (M e ) (Myr) (Myr) 



A Cas 


1.4 


0.88 


0.012 


0.36 


0.0034 


1.9 


7.3xl0" 4 


HR 333 




0.31 


0.040 






0.034 


7.8xl0- 3 


HR 506 


82 


0.24 


0.0063 


0.003 


0.13 


0.020 


1.1x10 s 


49 Cet 




0.15 


0.058 






0.0073 


9.3xl0~ 5 


7 Tri 




0.25 


0.019 






0.14 


3.0xl0~ 5 


t 3 Eri 




0.049 


0.012 






0.057 


3.1xl0~ 9 


HR 1082 




0.93 


0.061 






0.025 


1.3xl0 10 


HR 1570 




0.54 


0.013 






0.22 


2.0xl0 2 


C Lep 




0.031 


0.032 






0.011 


3.5xl0~ 12 


T T T ~\ — o 1 A o 

HD 53143 


o 
I 


0.014 


0.0013 






0.0036 


1.5xl0 -5 


HR 3314 




0.20 


0.018 






0.11 


3.9xl0- 6 


HD 95086 




0.37 


0.0071 






0.0081 


7.3xl0 3 


(3 UMa 




0.69 


0.026 






0.7 


2.7xl0~ 2 


(3 Leo 




0.12 


0.0049 






0.080 


3.8X10" 5 


HD 110058 






0.015 










A Boo 




0.16 


0.017 






0.085 


2.5xl0~ 5 


a CrB 




0.26 


0.031 






0.28 


3.4xl0- 7 


HD 139664 


? 


0.10 


0.0059 






0.058 


1.7xl0 _1 


HD 146897 


? 


0.13 


0.011 






0.00042 


1.4xl0 _1 


HR 6297 




0.17 


0.021 






0.027 


2.0xl0- 3 


HR 6486 




0.08 


0.0099 






0.026 


1.2xl0 -5 


HR 6532 


1.2 


0.48 


0.050 


0.17 


0.015 


0.2 


2.2xl0~ 4 



Table 5 — Continued 







PR Drag 


CPR Drag 


Collision 


Sublimation 


Name 


M ■ j(? 
wind 


tpR 


M PB 


tpR+wind 


Mpb 




f 






(Myr) 


(M e ) 


(Myr) 


(M e ) 


(Myr) 


(Myr) 


78 Her 




0.34 


0.0036 






0.29 


4.3xl0~ 7 


7 Oph 




0.20 


0.034 






0.052 


2.3xl0" 5 


HR 6670 


190 


0.081 


0.015 


0.00043 


0.69 


0.068 


1.9xl0 -3 


HD 181327 


? 


0.24 


0.99 


? 


? 


0.0013 


6.0xl0 4 


HD 191089 


? 


0.11 


0.65 


? 


? 


0.0014 


2. 7x10^ 


HR 8799 




0.027 


0.016 






0.014 


1.3xl0 -9 
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Table 6. X-ray and Stellar Wind Properties of FGKM Stars 



AT 

JName 


Onset 


Rub AT 


HR1 


7 1 


71 JT 




( ) 


(cts s J 




(erg s cm J 


( m q) 


107 PSC 


0.13 


027±0 003 


-0 92±0 04 


5.1xl0 3 


3.0 


HR 506 


0.46 


0.092 


-0 63±0 03 


1.7x 10 5 


290 


HD 16157 


0.10 


3 52±0 02 


11±0 01 






r 1 Eri 


0.04 


72±0 02 


-0 25±0 03 


l.OxlO 6 


? 


q/ For 


0.20 


2.8±0.2 


-0 08±0 05 


2.3xl0 6 


? 


'y Dor 


0.11 


19±0 04 


-0 09±0 18 


6.7xl0 5 


2100 


HR 1686 

-LJ--LIj 1UUU 


0.18 


020±0 007 


-1.0±0.4 


1.5xl0 4 


26 


HR 2483 




<0.05 




<1.5xl0 5 


<260 


HD 53143 


0.48 


19±0 01 


-0 37±0 05 


8.8x 10 5 


? 


HR 3220 


0.14 


35±0 03 


-0 18±0 09 


7.1xl0 5 


4100 


FI Vir 


0.05 


0.053 


-0.50±0.06 






HR 3862 


0.07 


0.002±0.0007 




<2.9xl0 3 


<2.1 


7) Crv 


0.02 


0.21±0.03 


-0.48±0.13 


3.4xl0 5 


1100 


HD 113766 




<0.05 




<7.6xl0 6 


? 


HR 5236 




<0.05 




<7.2xl0 5 


<2200 


a Boo 


0.11 


0.17±0.02 


-0.32±0.11 


2.5xl0 5 


650 


HD 139664 


0.11 


0.52±0.04 


-0.28±0.06 


1.2xl0 6 


? 


HR 5933 


0.35 


0.04±0.02 


-0.90±0.14 


1.2xl0 4 


14 


HD 146897 




<0.05 




<5.6xl0 6 


? 


HR 6670 


0.03 


0.10±0.02 


-0.50±0.17 


4.4xl0 5 


1700 


HD 181327 




<0.05 




<1.3xl0 6 


? 


HD 191089 


0.06 


0.02±0.02 


-0.30±0.22 


1.4xl0 6 


? 


HD 221354 




<0.05 




<3.7xl0 4 


<160 
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Table 7. Binary Properties 



Name 


Primary & Secondary 
Spectral Types 


Separation 
(") 


Separation 
(AU) 


Period 

(years) 


Notes f 


References 


A Cas 


B7V & B8V 


0.38 


41 


515 


VB 


2 


HR 1570 


AOV & ? 


40 


1500 






8 


HR 1686 


F6V & ? 


30 


630 






8 


k Lep 


B7V & B9V 


2.2 


370 






3,6 


32 Ori 


B5IV & B7V 


1.3 


120 


590 


VB 


1 


HR 2124 


A2V & ? 


0.1 


5 






5 


ip 5 Aur 


GOV & MOV 


31 


520 






8 


G Cen 


B3Vn & ? 


23 


2800 






8 


HD 113766 


F3/F5V & ? 


1.3 


170 






3 


HR 6532 


B9.5V & <A9V 






0.019 


SB 


4 


HR 7329 


AOVn & M7/8V 


4 


200 


1700 


CPM 


7 



tCPM = Common Proper Motion, SB = Spectroscopic Binary, VB = Visual Binary 



References. — (1) Abt & Cardona 1984; (2) Abt & Boonyarak 2004; (3) Fabricius & Makarov 
2000; (4) Guthrie et al. 1986; (5) Horsch et al. 2002; (6) Lindroos 1985; (7) Lowrance et al 
2000; (8) Simbad 
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Table 8. Emission Line and Gas Mass 3a Upper Limits 



50 K 100 K 



Name 


H 2 S(l) 




SI 




H 2 S(0) 




M H2 


M H2 




(erg s 1 cm 2 ) 


(erg s 1 cm 2 ) 


(erg s 1 cm 2 ) 


(M e ) 




49 Cet 


<8.6xl0" 


•15 


<1.0xl0" 


■14 


<7.4xl0~ 


•15 


<6000 


<80 


r Cet 


<2.7xl0" 


•13 


<3.2xl0" 


■14 


<1.9xl0" 


14 


<70 


<1 


r 3 Eri 


<5.0xl0~ 


■15 


<1.5xl0" 


-14 


<5.8xl0- 


■15 


<900 


<10 


a For 


<4.3xl0" 


•14 


<1.4xl0" 


-14 


<1.3xl0" 


14 


<600 


<7 


C Lep 


<2.1xl0" 


■14 


<1.8xl0" 


-14 


<1.0xl0" 


■14 


<1000 


<10 


HR 2124 


<2.3xl0~ 


•14 


<7.7xl0- 


-15 


<l.lxl0 _ 


14 


<5000 


<70 


5 Vel 


<5.0xl0- 


-14 


<1.2xl0" 


-14 


<4.2xl0" 


15 


<500 


<7 


(3 UMa 


<1.9xl0" 


■14 


<7.3xl0^ 


-15 


<3.8xl0~ 


■15 


<500 


<6 


j3 Leo 


<2.5xl0~ 


■14 


<6.1xl0" 


-14 


<4.4xl0- 


•14 


<1000 


<10 


i] Crv 


<2.3xl0~ 


14 


<1.2xl0" 


-14 


<8.9xl0- 


■15 


<600 


<8 


HD 113766 


<l.lxl0" 


■13 


<2.1xl0" 


-14 


<1.0xl0" 


14 


<40000 


<500 


7 Tra 


<2.7xl0~ 


14 


<1.4xl0" 


-14 


<9.3xl0- 


15 


<7000 


<80 


a CrB 


<3.0xl0- 


■14 


<1.3xl0" 


-14 


<6.5xl0~ 


15 


<800 


<10 


41 Ser 


<3.7xl0~ 


14 


<1.3xl0" 


-14 


<8.0xl0- 


■15 


<200 


<3 


HR 6486 


<2.6xl0" 


■14 


<1.0xl0" 


-14 


<9.0xl0- 


■15 


<1000 


<20 


7 Oph 


<1.3xl0" 


■14 


<1.9xl0" 


-14 


<6.4xl0- 


■15 


<1000 


<20 


HR 7012 


<2.4xl0" 


■14 


<l.lxl0" 


-14 


<1.3xl0" 


14 


<2400 


<30 
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Fig. 1. — Spectral Energy Distributions (SEDs) for all objects in our sample. TD1 fluxes 
(Thompson et al. 1978) are plotted as pentagons, General Catalogue of Photometric Data 
mean UBV or Johnson et al. (1966) fluxes are plotted as triangles, and 2MASS JHK fluxes 
(Cutri et al. 2003) are plotted as squares. Color-corrected IRAS, MIPS, and submillimeter 
photometry, where available, are shown with green, blue, and black error bars and upper 
limit symbols, respectively. Our Spitzer IRS spectra, as reported here, are shown in red. 
Overlaid are the best fit 1993 Kurucz and Nextgen models for the stellar atmospheres. 
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Fig. lb. - 
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Fig. lc. - 
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Fig. Id. - 
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Fig. le. - 
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X (|im) 



Fig. 2.— Photosphere-subtracted IRS spectra of (a) HR 3927, (b) rj Crv, (c) HD 113766, 
(d) HR 7012, and (e) r\ Tel. The solid magenta lines show the final models for each disk 
system. 
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Fig. 3a. — Photosphere subtracted IRS spectra with F v plotted as a function of wavelength. 
The minimum \ 2 fits for the single temperature black body and continuous disk models 
are overplotted in red and blue, respectively. IRAS photosphere subtracted fluxes, where 
available, are overplotted in green to highlight any discrepancies with IRS data. 
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Fig. 3b - 



-62 - 




Fig. 3c - 
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Fig. 4. — (a) The grain lifetimes are plotted as a function of distance around a B5V star. 
The Poynting- Robertson Drag/Stellar Wind drag lifetime is shown with a solid line; the 
sublimation lifetime is showed with a dotted line; and the collisional lifetime is shown with 
a dashed line, assuming M su j, mm = 0.001, 0.01, 0.1, and 1 M e (from top to bottom), (b) 
same as (a) for A5V star. No stellar wind drag is assumed, (c) same as (a) for a F5V star. 
The Stellar Wind drag lifetime is shown with a solid line, assuming that M wind = 100, 1000 
M Q (from top to bottom). 
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Fig. 5. — Histogram showing the distribution of inferred black body grain temperatures for 
objects whose IRS spectra are well-fit by a single temperature black body. 
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LOG Age (Myr) LOG Age (Myr) 

Fig. 6. — (a) Inferred fractional luminosity, L/^/L* plotted as a function of age. The solid 
lines show the Kenyon & Bromley (2005) model for the evolution of L/^/L* for dust at 3 - 
20 AU (left) and at 30 - 100 AU (right) around a main sequence A-type star. The dotted 
line shows the slope expected if the dust mass declines inversely with age as expected if 
collisions are the dominant grain destruction mechanism; HD 113766 and r\ Crv possess 
extremely high values of Lm/L* for their ages and define the fitting coefficient for the 1/t 
trend line. The dashed line shows the slope expected if Poynting-Robertson and corpuscular 
stellar wind drag are the dominant grain destruction mechanism; this trend line accounts for 
all of the stars in our sample except for rj Crv. (b) Black body grain distance plotted as a 
function of stellar age. Stars with proposed multiple debris belts are shown in color with a 
line connecting multiple points. HR 3927, r\ Crv, HD 113766, HR 7012, and r] Tel are shown 
in red, green, blue, cyan, and magenta, respectively. In both plots, B-type stars are shown 
with filled triangles; A-type stars are shown with filled squares; F-type stars are shown with 
filled circles; the K-type star is shown with a cross. The addition of submillimeter data to 
infer L IR /L* and D for 77 Crv and r Cet (shown as open circles) does not appear to elucidate 
either the relation between fractional infrared luminosity and age or grain distance and age. 
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Fig. 7. — Inferred fractional luminosity, Ljr/L*, and grain temperature, T gr , plotted as 
a function of binary separation (for all binaries in our study), measured stellar rotational 
velocity, v sini, and stellar metallicity, [Fe/H]. For binary systems to the left of the dotted 
line, the LL (and sometimes the SL) slit contain both the primary and secondary; for objects 
to the right of the dotted vertical line, the SL and LL slits contain only the primary star. In 
plots of grain properties as a function of v sin 2 and [Fe/H], B-type stars are shown as crosses, 
A-type stars are shown as squares, F-type stars are shown as circles, and K-type stars are 
shown as triangles. Stars with spectral type earlier than F5V are shown as open symbols 
in v sin % plots. Stars with spectral type earlier than F0V are shown with open symbols in 
[Fe/H] plots. 



